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Abstract 


A  review  of  theoretical  developments  in  predicting  the  buckling  response  of  cylinders  subject 
to  impulsive  loads  is  presented.  Most  of  this  theory  dealc  with  axisvmmetric.  radial  impulses 
on  cylinders.  The  development  of  solutions  for  the  critical  modes  and  loading  magnitudes 
which  produce  excessive  growth  of  displacements  are  reviewed.  Existing  theories  cover  the 
specific  cases  of  either  entirely  elastic  or  entirely  plastic  material  behaviour  for  infinite  length 
and  short  cylindrical  shells.  The  resultant  theories  are  applied  to  various  shell  geometries  to 
investigate  influencial  parameters.  A  review  of  numerical  finite  element  and  finite  difference 
studies  which  investigate  dynamic  pulse  buckling  is  also  given.  The  requirements  to  examine 
dynamic  buckling  of  more  complex  structures  such  as  ring  stiffened  cylinders  are  discussed. 


Resume 

Un  examen  des  perfectionnements  theoriques  pour  la  prediction  de  la  reaction  au  flambage 
des  cylindres  soumis  a  des  impulsions  de  charge  est  presente.  La  plus  grande  partm  de 
cette  theorie  a  pone  sur  les  impulsions  radiales,  axisymetriques  exerc6es  sur  les  cylindres. 
On  y  examine  aussi  les  solutions  mises  au  point  pour  les  amplitudes  de  charges  et  les 
modes  critiques  qui  produisent  des  deplacements  sans  limites.  Les  theories  existantes 
portent  sur  des  cas  precis  de  comportement  des  materiaux  entierement  elasriques  ou 
entierement  plastiques  pour  des  enveloppes  cylindriques  courtes  et  des  longueur  infinie. 

Les  theories  qui  en  decoulent  sont  appliquees  a  differentes  geometries  d'enveloppes  afin 
d’etudier  les  parametres  d'influence.  Un  examen  des  etude  d'elements  finis  numeriques  et 
decarts  finis  relatives  a  l'etude  du  flambage  par  impulsion  dynamique  est  aussi  presente. 
On  y  traite  des  exigences  pour  l'etude  du  flambage  dynamique  de  structures  plus  complexes 
comme  des  cylindres  renforces  de  frettes. 


Contents 

Abstract 

Table  of  Contents 
List  of  Figures 
Notation 

1  INTRODUCTION 


Aooasalon  Pc 
KTIS  GRAM* 
MIC  TAB 

Ufcaanounoed 
I  Juatlfioatlo 


By _ _ _ 

Distribution/ 
Availability  Codes 
j Avail  asd/#r 
Dlst  I  Special 

L 


f- 


ii 

iii 

iv 

v 

1 


2  THEORETICAL  REVIEW  3 


2.1  Plastic  Flow  Pulse  Buckling .  5 

2.1.1  Derivation  of  the  Equation  of  Motion  .  5 

2.1.2  Particular  Solution .  7 

2.1.3  Perturbed  Velocity  Solution .  8 

2.1.4  Complete  Solution .  9 

2.1.5  Discussion  of  Plastic  Flow  Buckling  .  10 

2.2  Elastic  Pulse  Buckling .  14 

2.2.1  Derivation  of  the  Equation  of  Motion  . 15 

2.2.2  Impulse  Velocities  for  Vibration  and  No  Buckling .  17 

2.2.3  Pulse  Buckling  from  Large  Initial  Velocities . 22 

2.3  Elastic  and  Plastic  Pulse  Buckling . 27 

2.4  Pulse  Buckling  for  Finite  Length  Cylinders  . 31 

2.4.1  Plastic  Flow  Buckling  of  Finite  Length  Cylinders . 31 

2.4.2  Elastic  Pulse  Buckling  of  Finite  Length  Shells . 42 

2.5  Effect  of  Pulse  Duration  on  Response . 46 

2.6  Effect  of  Spacial  Pulse  Shape . 51 

2.7  Strain  Rate  Reversal . 52 

2.8  Material  Property  Effects . 54 


3  Approximate  Formulae  for  Critical  Impulse  and  Critical  Modes  of  Dynamic 

Pulse  Buckling  57 

4  Review  of  Numerical  Solutions  to  Dynamic  Pulse  Buckling  64 

5  Conclusions  68 


APPENDIX  A:  Derivation  of  the  Shell  Curvature  Expression 


iii 


69 


APPENDIX  B:  Derivation  of  the  Mathieu  Stability  Equation  71 

References  73 

List  of  Figures 

1  Strain  Differential  from  Strain  Hardening  (from  reference  [9]) .  6 

2  Shell  Force  Components .  7 

3  Modal  Growth  for  Model  2  .  12 

4  Modal  Growth  for  Model  9  .  13 

5  Mathieu  Stability  Diagram  (from  reference  [10]) .  18 

6  Response  Showing  Transfer  of  Energy  From  Hoop  to  Flexural  Modes  vfrom  ref¬ 
erence  [4] ) .  ■  r‘ 

7  Elastic  Motion  for  Various  Stability  Parameters,  p  (from  reference  [4]) . 22 

8  Mathieu  Stability  Diagram  for  Large  Initial  Velocities  Causing  Buckling  (from 

reference  [13]) .  23 

9  Modal  Amplification  for  Elastic  Dynamic  Buckling  of  Thin  Walled  Cylinder  (from 

reference  [13]) .  24 

10  Hoop  and  Flexur.  1  Yield  and  Dynamic  Buckling  Limits  for  Elastic  Response  to 

Pulse  Loads  (from  reference  [4]) .  25 

11  Strain  Vectors  on  Yield  Ellipse  (from  reference  [20] )  35 

12  Displacement  and  Velocity  Amplifications  as  a  Function  of  Harmonic  Number 

for  A  Short  Cylinder  (from  reference  [20]) .  37 

13  Yield  Ellipse  for  Variable  Length  Shells  (from  reference  [21]) . 38 

14  Mathieu  Stability  Curve  for  Finite  Length  Shell  (from  reference  [23]) . 43 

15  Peak  Pressure  versus  Impulse  Curves  (from  reference  [19]) . 47 

16  Peak  Pressure  versus  Impulse  Curves  (from  reference  [19]) . 48 

17  Critical  Pulse  Parameters  for  Axisymmetrically  Loaded  Cylinders  ( from  reference 

[’25]) .  50 

18  Comparison  of  Finite  Element  Results  Including  Strain  Rate  Reversal  with  Sim¬ 
pler  Plastic  Flow  Results  (from  reference  [4]) .  53 

19  Material  Property  Curve  (from  reference  [4] )  56 

20  Critical  Impulse  Versus  a/h  Ratio  (from  reference  [4])  . 60 

21  Static  and  Dynamic  Response  of  a  Perfect  and  Imperfect  Sphere  (from  reference 

[29]) .  65 


IV 


Notation 

a 

MrKMT) 


do  -  u'o 


A 


ncr 


An(0-Bn(0 

Bn 

c 


Ch 


shell  radius  to  mid  thickness 
asymmetric  amplitudes  of  motion 
radial  hoop  mode  amplitude 

displacement  amplification  from  autoparametric  excitation  of  flexural 
mode 

amplification  functions  for  displacement  and  velocity  imperfections 
=  v0an/pn 

speed  of  sound  in  material  =  \T^f  *  constant  in  equation  (159) 

=  \rBhTp 


C(q).  f.q.qo.c.n  coefficients  defined  in  equation  (161) 
Ci .  Co.  C'3,  C4.  C5.  C$  constants  defined  in  reference  17 
d A  differential  arc  length  of  shell 


D 

Eh 

Ei 

E 

E, 

fn- 

F 

h 

ho 

/ 


—  t’O '  /  Pn 

12  E 
(l-^2) 

strain  hardening  tangent  modulus 
equivalent  Young's  Modulus  = 
Young's  Modulus 
tangent  modulus 

generalized  variables,  dependent  on  « 

stress  function 

shell  thickness 

shell  thickness 

moment  of  inertia 


v 


1a 

hi 

Ic 

k 

A'i 

A'2 

a-3 

K 

L 

l 

Mi 

m 

n 

71  cr 

y, 

Pn 

P 

P 

Po(  T  ) 
Pn(r) 

Pa 

Po 

Pn 

P-D 


asymptotic  peak  impulse 
critical  impulse  from  Dirac  delta  function 
peak  impulse  amplitude 
variable  length  parameter 
=  2(1  -  k  +  k2) 

(3A'1/2)1/2 
(2  —  k)/ K'i 

material  curve  parameter  defined  as 
length  of  shell 

nondimensional  length  =  L/a 

resultant  moment  in  i  direction 

mass  of  shell  material  =  ph 

circumferential  harmonic  wave  number 

circumferential  harmonic  of  greatest  unbounded  growth 

resultant  membrane  force  in  i  coordinate  direction 

=  (n2  —  l)(n2  -  s2  ) 

stability  parameter  =  V'0/co 

radial  pressure 

axisymmetric  pressure  amplitude 
asymmetric  pressure  amplitude 
asymptotic  peak  pressure 
minimum  pn'se  amplitude  to  cause  buckling 

_  4q2~»V’q 

3(2  —  k)a 

constants  defining  material  curves  that  are  strain  rate  dependent 


VI 


Q 

Qn 

T 

Rl 

R 

S 

s2 

Sn 

t 

T 

Tt 

Te 

h 

u 

V 

~  Ji 

«c 


U 

V'o 


»'0 


r 


shear  force 

parameter  of  Bessel  equation,  defined  in  text 
deformed  shell  radius  =  a  -  iv0 
parameter  of  Bessel  equation,  defined  in  text 
radius  of  deformed  motion,  =  a(  1  —  it) 
circumferential  membrane  force 

nondimensional  shell  parameter=  Sa2/Ehl  =  \2aa2 / E^h2  for  plastic 
flow  buckling. =  p/a  for  elastic  buckling 

parameter  of  Bessel  equation,  defined  in  text 

time 

kinetic  energy  or  pulse  duration  time 

pulse  time  defining  plastic  flow  solution 

pulse  time  defining  elastic,  quasi-static  solution 

time  spent  in  elastic  motion 

nondimensional  radial  displacement  =  it /a 

strain  energy  of  shell 

si  ra i n  p*- -rtrv  of  Iiot  mode 

strain  energy  of  flexual  mode 

complimentary  solution  of  flexural  mode 

modal  amplitudes  for  axial  harmonic,  in  and  circumferential  harmonic, 
n 

particular  solution  to  hoop  mode 
initial  radial  velocity 

nondimensional  initial  velocity  =  \f\2p/ Eha\'u[ h 
tangential  displacement 

vii 


1C 

U', 

U'n 

u'y 

u'J 

Z 


Qn  .  Jn 

Q 2 

3 


( 

£mi 
(  f 


(. 

Tr. 

1 


1  r 


radial  displacement 
initial  shape  imperfections 
natural  frequency  of  vibration 
radial  motion  in  elastic  regi m<> 
final  radial  motion 
radial  inertia  force 
coordinate  through  shell  thickness 
perturbations  of  the  initial  velocity 
hi 

~  1  2a  ‘ 

critical  damping  ratio,  or  =  \J Eh/ 1' 
amplitudes  of  initial  imperfections 
nondimensional  initial  shape  imperfection  =  j-<*n 

strain 

strain  of  shell  midsurface 

midsurface  strain  of  initial  imperfections 

flexual  strain 

strain  at  yield 

generalized  strain  rate 

—  pXVl  —  A2 

material  viscosity  constant 
perturbation  amplitude  of  sine  function 
curvature  for  initial  shape  imperfections 


viii 


shell  curvature 


A 


Pn 

V 


o 

V 


Up 

p 

(T,,, 


Vf{ 

af 


a 

a! 


vo 

To 


‘  f 

T±nn 

8 


wave  parameter  =  n/s 

parameter  of  Mathieu  instability  equation 

Poisson's  Ratio 

parameter  of  Mathieu  instability  equation 
angle  of  curvature  after  deformation  ^  ^ 

nondimensional  tangential  displacement  =  v/a.  also  change  in  shell  angle 
of  curvature.  8  —  0 

deformed  radius  of  curvature 

shell  material  density 

constant  plastic  flow  stress 

maximum  stress  in  hoop  mode 

maximum  stress  in  flexural  mode 

yield  stress 

generalized  stress 

deviatoric  stress  of  component  i 

current  midsurface  stress 

generalized  stress  at  infinitesimally  small  strain  rate 
termination  time  of  elastic  motion 
nondimensional  time,  defined  in  various  sect  ions 
time  to  terminal  radial  motion 
time  at  which  strain  rate  reversal  occurs 

circumferential  coordinate,  also  initial  angle  of  curvature  of  the  shell 


1  INTRODUCTION 


Buckling  under  static  loading  is  traditionally  investigated  in  the  design  of  shell  structures. 
Buckling  instability  is  also  a  possible  mode  of  failure  in  structures  subject  to  dynamic  loads. 
Dynamic  buckling  of  impulsively  loaded  shell  structures  occurs  through  unstable  growth  of 
displacements  during  motion.  Dynamic  buckling  has  been  considered  in  the  design  of  such 
structures  as  aircraft,  automobiles,  trains,  pipelines  and  reactor  pressure  vessels  to  improve 
their  resistance  to  accidental  impact  or  explosion  loading.  A  dynamic  buckling  mode  of  failure 
provides  a  high  energy-absorbing  mechanism  which  decreases  the  transmitted  shock  to  internal 
components  or  personnel.  Improvement  in  design  to  resist  shock  loads  from  hostile  weapons 
or  operational  conditions  has  been  the  goal  of  considering  dynamic  buckling  in  the  design  of 
military  structures  such  as  missiles  and  rockets.  This  study  presents  a  review  of  theoretical 
developments  in  dynamic  pulse  buckling  to  delineate  the  physical  concepts  of  the  process.  This 
review  serves  as  a  basis  for  further  work  in  developing  an  understanding  of  dynamic  buckling 
in  the  context  of  submarine  pressure  hull  response  to  underwater  shock  loading. 

Studies  of  submarine  response  to  shock  loading  have  concentrated  on  predicting  the  stress 
level  and  pattern  resulting  from  complex  pulse-structure  interaction  [1.2]  and  have  generally  not 
considered  dynamic  buckling  failure  which  will  probably  occur  after  the  material  has  reached 
its  yield  limit.  As  in  static  load  studies,  critical  combinations  of  structural  parameters  and 
load  will  cause  instability,  and  as  in  static  analysis,  detection  of  these  dynamic  stability  limits 
requires  special  formulae  and/or  numerical  methods.  Previous  work  in  dynamic  buckling  has 
been  mainly  for  axisymmetric.  pulse  loading  of  unstiffened  cylinders.  The  submarine  problem  is 
more  <~omplex  as  it  requires  consideration  of  asymmetric  pulse  loading  of  ring  stiffened  cylinders. 

The  term  "dynamic  buckling'  has  been  used  to  describe  two  different  types  of  structural 
behaviour  mechanisms.  The  first  type  is  buckling  of  a  structure  subject  to  periodic  loading 
functions  which  result  in  a  resonance  with  a  buckling  mode.  This  is  termed  "parametric  buck¬ 
ling-.  as  the  loading  function  is  a  parameter  of  the  displacements  in  the  differential  equation 
which  describes  the  motion  [3].  The  load  intensity  to  cause  collapse  can  be  lower  than  the 
static  buckling  load  for  parametric  buckling  behaviour.  The  second  type  of  dynamic  buckling 
is  usually  termed  "pulse  buckling'  resulting  from  a  transient  loading  function  of  a  single  pulse 
form.  Pulse  loading  intensities  needed  to  cause  buckling  are  larger  than  the  static  buckling 
load.  This  review  is -concerned  primarily  with  the  latter,  ‘pulse  buckling'  definiton  of  dynamic 
buckling,  as  it  corresponds  to  buckling  caused  by  shock  loads. 

An  extensive  review  of  the  theory  and  analytical  solutions  to  the  pulse  buckling  problem  is 
given  in  Section  2  for  the  purpose  of  developing  an  understanding  of  the  physics  of  dynamic 
pulse  buckling.  The  Stanford  Research  Institute  (SRI)  has  devoted  effort  to  this  topic  for  the 
last  three  decades  and  most  relevant  work  comes  from  this  source.  A  collection  of  SRI  work 
has  been  compiled  in  a  comprehensive  manuscript  [4].  Other  reviews  of  dynamic  pulse  buckling 
are  given  in  References  [5,6.7],  As  is  the  case  for  theoretical  studies  of  most  complex  problems. 
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analytical  solutions  have  been  derived  only  for  relatively  simple  geometries  and  loading  func¬ 
tions.  Asymmetric  loading,  ring  stiffeners  and  other  con  •  lexities  make  analytical  solutions  for 
pulse  buckling  of  submarine  pressure  hulls  unlikely.  Numerical  modelling  by  finite  element  or 
finite  difference  methods  offers  a  means  of  solving  these  more  complex  problems.  In  any  type 
of  analysis,  one  must  kn«  w  before  hand  what  structural  behaviou  uo  expect,  and  formulate 
a  model  accordingly.  The  physics  of  dynamic  buckling  is  complex,  in  that  one  may  not  know 
if  the  response  will  be  elastic  or  plastic,  of  short  duration  or  long,  a.  d  indeed,  one  may  also 
not  be  able  to  define  the  mode  or  point  of  failure.  Therefore  an  understanding  of  the  theory 
of  dynamic  buckling  for  simple  cases  is  esajntial  before  undertaking  more  complex  numerical 
investigations. 

Section  3  of  this  report  gives  approximate  formulae  for  the  determination  of  critical  buckling 
modes  and  critical  impulse  loads  for  dynamic  buckling  of  simplified  cases.  These  have  resulted 
from  assumptions  made  to  the  various  theories.  The  effects  of  several  shell  geometry  and 
material  parameters  on  the  critical  modes  and  loads  are  investigated  through  those  formulae. 

A  few  studies  have  investigated  response  of  shells  to  impulsive  loading  by  finite  element  or 
finite  difference  methods.  A  review  of  some  of  this  work  is  given  in  Section  4. 
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2  THEORETICAL  REVIEW 


In  static,  linearly  elastic,  bifurcation  buckling,  the  critical  load  is  a  well  defined  point 
depending  only  on  the  structural  geometry,  material  properties  and  loading  distribution.  In 
dynamic  pulse  buckling,  the  solution  is  for  excessive  growth  of  a  particular  mooe  or  modes, 
rather  than  for  specific  critical  points.  The  intensity  and  duration  of  the  load  affect  t  h*>  buckline 
mode  as  well  as  the  total  response.  The  response  also  depends  on  the  time  spent  in  the  elastic 
and  plastic  phases  of  the  motion. 

This  section  reviews  theories  for  simplified  physical  cases  of  plastic  flow  pulse  buckling  and 
elastic  pulse  buckling.  All  theories  start  with  the  derivation  of  the  equations  of  motion  for  the 
physical  problem  of  concern.  They  then  proceed  to  the  solution  for  the  unbuckled,  axisym- 
metric.  radial  hoop  motion.  The  complimentary  solution  for  the  flexural  buckling  motion  is 
then  considered  for  initial  displacement  and  velocity  perturbations  with  harmonic  circumferen¬ 
tial  variation.  Nonlinear  ordinary  differential  equations  result,  requiring  numerical  integration 
for  complete  solution,  unless  linearizing  assumptions  can  be  justified.  Dynamic  instability  is 
determined  by  investigating  the  boundedness  of  the  perturbed  motion;  this  is  the  Liapunov 
definition  of  instability  [8].  In  all  cases  investigated,  the  nonlinear  differential  equations  am 
linearized  and  examined  qualitatively  to  determine  their  stability  via  exponential  growth  of 
certain  critical  modes. 

Solutions  for  plastic  behaviour  (thick  shells  with  small  radius  to  thickness,  a/h.  ratios)  and 
elastic  behaviour  (thin  shells  with  large  a/h  ratios)  have  been  derived  separately  [9.10]  and 
are  reviewed  here  to  illustrate  the  governing  physics  of  each  case.  Basic  assumptions  of  either 
entirely  plastic  or  entirely  elastic  behaviour  are  needed  to  permit  solution.  The  critical  radius 
to  thickness  ratio  for  transition  from  predominantly  elastic  to  predominantly  plastic  response 
is  in  the  range  of  100  to  200  for  most  engineering  materials.  The  transition  is  not  distinct. 
Completely  plastic  behaviour  can  be  assumed  for  a/h  ratios  below  40  and  completely  elastic- 
behaviour  can  be  assumed  for  a/h  ratios  greater  than  400.  Shells  of  intermediate  a/h  ratios  have 
some  combination  of  plastic  and  elastic  behaviour.  The  magnitude  of  the  pulse  also  affects  the 
ultimate  response.  Solutions  for  intermediate  a/h  ratio  shells  require  consideration  of  significant 
contributions  from  both  elastic  and  plastic  response,  and  no  known  analytical  solutions  have 
been  derived  for  this  case. 

This  section  presents  discussions  of  the  individual  plastic  and  elastic  theories  followed  by  a 
discussion  of  the  work  of  Stuiver  [1 1]  who  presents  a  common  method  of  deriving  both  the  elastic 
and  plastic  equations  and  shows  that  the  plastic  response  has  some  dependency  on  the  elastic 
motion.  Other  studies  which  investigate  the  effects  of  shell  length,  pulse  duration  and  spacial 
shape,  varying  material  properties  and  strain  rate  reversal  on  the  shell  motion  and  stability  are 
then  reviewed.  In  order  to  facilitate  the  use  of  the  various  theories  in  design,  most  authors  have 
applied  simplifying  assumptions  to  derive  approximate  formulae  for  critical  dynamic  buckling 
modes  and  impulses.  These  are  presented  in  Section  .3  along  with  sample  calculations  using 


these  formulae  for  a  wide  range  of  shell  parameters. 
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2.1  Plastic  Flow  Pulse  Buckling 

Abrahamson  and  Goodier  [9],  produced  the  first  satisfactory  solution  to  the  dynamic  pulse 
buckling  of  shells  in  1962.  The  solution  was  for  completely  plastic  response  of  an  infinite  length 
cylinder  to  a  perturbed  axisymmetrir  velocity  pulse.  It  was  limited  to  cylinders  of  smaller  a/h 
ratios  where  elastic  behaviour  can  be  ignored.  This  is  the  simpler  of  the  theories  for  dynamic 
pulse  buckling  and  illustrates  the  physics  of  problem  very  well.  This  theoretical  work  was 
verified  by  considerable  experimental  investigation. 

The  basic  assumption  of  this  theory  is  that  motion  occurs  entirely  in  the  plastic  material 
regime  and  that  the  entire  shell  is  in  a  state  of  continuously  increasing  compressive  plane  strain 
with  no  strain  rate  reversal.  Perturbed  inextensional  flexural  modes  are  superimposed  on  the 
axisymmetric.  constant  compressive  plastic  flow  stress  state.  The  resistance  of  the  shell  to  the 
motion  comes  from  moments  produced  by  a  stress  differential  across  the  shell  thickness  as  a 
result  of  the  strain  hardening  or  tangent  modulus.  Thus,  this  theory  has  also  been  called  the 
"tangent  modulus’  theory  as  its  existence  is  dependent  upon  the  material  having  a  non-zero 
plastic  tangent  modulus. 

2.1.1  Derivation  of  the  Equation  of  Motion 

Figure  1  illustrates  the  stress  differential  occuring  in  a  shell  perturbation  (local  disruption  in 
uniform  radial  motion)  as  a  result  of  the  strain  hardening  modulus.  Eh-  This  stress  differential 
results  in  a  bending  moment: 


M  —  Eh  I  k  (1) 

where  I  =  pj h 3  for  a  unit  width  of  shell  and  K  is  the  curvature  of  the  shell  segment,  the 
derivation  of  which  is  given  in  Appendix  A.  and  is  defined  as: 


K 


1  02  w 
a 2  d62 


+  tr) 


(2) 


The  assumption  is  made  that  the  solution  is  for  a  ring  or  an  infinite  length  cylinder;  no 
end  effects  or  biaxial  stress  state  are  included  in  the  curvature  derivation.  The  radius  and  the 
thickness  of  the  shell  are  conservatively  assumed  to  remain  at  their  initial  value  throughout  the 
motion.  In  reality,  the  a/h  ratio  of  the  shell  will  decrease  with  the  motion  making  the  shell  more 
resistive  to  buckling.  The  mean  circumferential  hoop  stress.  am-  is  assummed  to  be  constant 
around  the  circumference  and  with  time.  In  reality,  am  will  vary  as  the  radial  deflection.  u.'(<p). 
varies  and  more  important,  it  will  also  increase  due  to  strain  hardening.  The  strain  hardening 
modulus  is  assumed  to  be  constant  and  crm  is  taken  as  the  average  stress  in  the  post  yield, 
plastic  flow  region.  This  assumption  affects  the  solution  by  overestimating  the  time  required 


Figure  1:  Strain  Differential  from  Strain  Hardening  (from  reference  [9]) 


for  the  hoop  motion  to  absorb  the  kinetic  energy  of  the  impulse.  The  average  compressive  hoop 
force  in  the  shell  is: 

S  h  (y  rri 


From  Figure  2.  equilibrium  of  forces  results  in: 

J2  M  =  (M  +  dM  )dx  -  Mdx  -  Q(dX)dx  =  0 

which  reduces  to: 


Q 


dM 
OX  ' 


and: 


which  reduces  to: 


T  Fz  =  ~{Q  +  dQ)dx  +  Qdx  +  ZdXdx  —  Sd<t>dx  =  0 


dQ  do  d2u • 

0A+SM  =  ~"'W 


3) 


(4) 


(6) 


(") 


where  m  is  the  mass/unit  circumference  and  where  If  =  -  in  which  -  is  the  deformed  radius 
of  curvature  defined  by: 

11  11  d2w 

-  =  -  +  K=  -  +  (8) 
pa  a  a~  dO- 
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Figure  2:  Shell  Force  Components 

Substituting  equations  (2)  and  (1)  into  equation  (5)  and  equations  (8)  and(5)  into  equation 
gives  the  partial  differential  equation  of  motion  for  the  shell: 

Ehld^u'  d2w  .1  1  d2w  .  d2u 

+  inrr>  +  SI-  +  -Aw  +  ”')]  =  -m 


q4  'as4  as1 


Reference  [9]  rewrites  equation  (9)  in  dimensionless  form  using: 

1 


ir  Ehl ,  1  Eh  ht  2 

u  =  — . r  =  \  - -t  -  —7=\ - =•  and  s  - 

a  V  m(,A  vl2  V  p  a2 


dt 2  ' 

Sa 2  12crtno2 


to  give  equation  (9)  as: 


d4  u  2 

w  +  {l  +  s 


d2u  2  d2u 

]d(p  +  s  u+ 


Ehl 


—  s 


EhU 


(9) 
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2.1.2  Particular  Solution 

The  particular  solution  to  equation  ( 1 1 )  is  for  the  unperturbed,  axisymmetric.  radial  motion 
of  the  perfect  ring  (ie.  u  is  independent  of  6).  Equation  (11)  is  reduced  to: 


d2u 

dr2 


+  s*u  =  -s 


I 


(12) 


wit li  initial  conditions  for  an  initial  velocity  pulse: 


dw  du 

u'(o)  -  ulo)  =  0.  and  (  — —  )(=o  =  1  o-  which  becomes  (  — )r=0=  r0 

at  or 


and  solution: 


up(r)  =  —  1  +  cossr  + 


Vo  . 

—  sin  sr. 
£ 


(1-1) 


The  duration  of  t lie  inward  hoop  motion  is  established  from  the  particular  solution  by  deter¬ 
mining  the  time  at  which  the  radial  velocity  vanishes.  Setting  the  first  derivitive  of  equation  (11) 
equal  to  zero  gives: 

1  .  r*  ,  i  -  • 

Tf  =  -  arc  tan  —  (lot 

s  s 


2.1.3  Perturbed  Velocity  Solution 

The  complimentary  solution  of  equation  (11)  is  derived  assuming  a  harmonic  circumferen¬ 
tial  perturbation  in  the  initial  velocity  profile  of  the  shell: 

Du  °° 

(  — )r  =  U  =  i'o [  1  +  ^(Qr;  COS  uO  +  ,3„  sill  1x6))  (10) 

°T  n  =  2 

where  an  and  3n  are  parameters  dependent  on  the  degree  of  velocity  perturbation.  This  results 
in  a  solution  of  the  form: 


u(t)  =  ]T[/n(r)cos7!0  +  ,(7n ( r ) si n  n9)  (17) 

n—2 

which  after  substitution  into  the  homogeneous  form  of  equation  (11)  gives: 

n4/ri(r)  cos  n6  +  <7n ( t ) sin  ixd  +  (1  +  s2)(-n2  fn(r)  cos  nO  -  n2gn(r ) sin  n9) 

+s2[Jn{r)  cos  n6  +  gn{T)s\nnO)  +  /n(r )  cos  nO  +  £,(r)  sin  =  0  (IS) 

yielding  the  ordinary  differential  equation: 

In  +  [n4  -  TX2(  1  +  .s2)  -I-  .s2]/n  =  o  (19) 

and  a  similar  equation  for  gn.  The  coefficient 

[n4  -  u2(l  +  ,s2)  +  .s2]  =  (u2  -  1  )(n2  -  s2)  =  p2  (20) 

determines  the  stability  of  motion.  For  n  <  s.  f„  and  <)n  are  hyperbolic,  unstable  functions  and 
for  a  >  s,  /„  and  gn  are  circular,  stable  functions.  In  other  words,  only  circumferential  modes. 


n,  which  are  less  than  the  shell  parameter,  s,  will  grow  exponentially  if  the  load  intensity  is 
great  enough.  For  example,  a  shell  of  parameters  a/h=30,  Eh  =  1  O'’  psi  and  rrm  =  75.000  psi  has 
a  value  5=29  which  allows  a  large  number  of  harmonics,  n,  to  grnw  unbounded.  One  of  these 
unbounded  harmonics  will  predominate. 

Taking  the  value  of  s  to  be  its  nearest  greater  integer  value,  the  complimentary  solution  of 
fn  and  gn  is: 

oo 

fn(r)=  y^(An  cosh pnT  +  Dn  si nh  p,,  r)  cos  nS  +  ^  ( Ari  cos pnT  +  Bn  sm  pnr  )co.-,n8  (21) 

n=2  n=5-f-l 

where  p2  =  (n2  -  l)(s2  -  n2)  for  n  <  s,  or.  =  (n2  -  1 )( r?2  -  .s2)  for  n  >  s  +  1. 


2.1.4  Complete  Solution 

Adding  equations  (21)  and  (14)  gives  the  complete  solution  as: 


u(t.  6)  =  -1  +  cos  st  +  —  sin  st  +  [(.4n  cosh  pnr  +  Bn  sinh  pnr)  cos  nd 


n=2 

OO 


+  (Cn  cosh  pnr  +  Dn  sinh  pnr)  sin  n9)  +  ^  [(A„  cos pnr  +  Bn  sin  pnr)  cos  n8 

«=j+ j 

+(C„  cos pnr  -f  Dn  sin  pnr)  sin  n0] 


(22) 


An  and  Cn  are  zero  from  the  initial  conditions  (13)  and  Bn  and  Dn  are  determined  from 
the  initial  velocity  perturbation,  equation  (1G),  to  be:  Bn  =  and  Dn  -  giving  the 
final  solution  as: 


u(r,  6)  =  - 1  +  cos  st  -f  ^  sin  st  +  r0  P)n=2(Q,'>  cos  +  4r»  sin  T>^)~  sinh  pnr 

+  t’o  En=s+i(Qn  cos n(9  +  !3n  sin  sil1  Pn~ 


(2.3) 


The  buckling  instability  occurs  as  a  result  of  excessive  growth  of  the  —  sinh  pnr  term.  One 


value  of  pn  will  dominate  the  response.  Taking  the  derivitive  of -V  sinh2  pnr  gives  the  limit  point 

rn 


as: 


l)2.  or.  n2r  -  i(.e2  +  J 


(24) 


For  the  shell  dimensions  used  in  the  above  example,  ( pn)cr  is  420  giving  an  nCT  of  21  as  the 
critical  mode  for  dynamic  pulse  buckling.  The  function  j^sinhpn7-  should  be  investigated  for 
a  range  of  n  in  the  vicinity  of  ncr,  and  for  a  duration  tj  defined  by  equation  (15)  to  determine 
the  instability  characteristics. 
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Since  this  theory  has  used  the  assumption  of  monotonicall y  increasing  compressive  strain, 
the  limit  point  of  strain  rate  reversal  has  to  he  investigated.  The  total  compressive  strains 
consisting  of  the  hoop  strain  and  flexural  strain  for  the  inside  and  outside  surfaces,  respectively, 
are: 


e(r.e)  =  ii  +  i[-  + 


h  02w,  ,  h  u- 


h  d2  w 


2a  a  2aL  (Jd‘  '  ‘  2a"  a  2a 2  dO2 

These  equations.  upon  substitution  of  equation  (23)  have  to  be  investigated  at  various  times 
throughout  the  motion  to  determine  the  onset  of  strain  rate  reversal.  The  occurrence  of  strain 
rate  reversal  is  an  indication  that  instabilitv  is  imminent. 


2.1.5  Discussion  of  Plastic  Flow  Buckling 

A  number  of  assumptions  have  been  made  in  Abrahamson  and  Goodier’s  plastic  flow- 
buckling  theory.  They  are: 

1.  The  solution  is  for  plane  strain,  and  therefore  an  infinite  cylinder  or  ring. 

2.  The  loading  function  is  a  perturbed,  axi symmetric,  radial  velocity  impulse  at  the  shell 
surface. 

3.  The  solution  ignores  the  elastic  behaviour  completely,  and  therefore  is  applicable  to  shells 
with  small  a/h  ratios  in  which  significant  plastic  flow  will  occur. 

■1.  The  radial  hoop  motion  is  assumed  to  occur  with  a  constant  stress  value,  and  with  const  ant 
initial  values  of  shell  thickness  and  radius.  The  tangent  modulus  is  also  assumed  to  be 
constant  during  the  motion. 

5.  The  stress  state  is  one  of  continuously  increasing  compression  and  therefore,  the  theory 
is  only  applicable  up  to  the  point  of  strain  rate  reversal. 

All  of  these  assumptions  have  been  addressed  in  studies  subsequent  to  that  of  Abrahamson  and 
Goodier  [9],  which  are  discussed  in  later  sections. 

Defining  the  velocity  pulse  parameters,  r,,on  and  0„  is  the  main  difficulty  in  using  this 
method  to  define  a  structure's  response  to  a  shock  load.  Reference  [9]  investigated  several 
velocity  pulse  profiles  and  although  the  buckled  shapes  differed,  the  results  only  varied  by  5 
percent. 

For  given  shell  dimensions,  thin  theory  may  be  used  to  establish: 

•  the  expected  dominant  mode  shape  of  plastic  pulse  buckling  via  equation  (24). 

t  an  approximation  to  the  duration  of  the  response  time  via  equation  (15). and 

•  the  motion  of  the  shell,  given  the  velocity  pulse  parameters  via  equation  (231. 
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To  determine  the  critical  impulse  for  dynamic  buckling,  one  must  establish  a  limit  point  at 
which  buckling  is  assumed  to  occur.  This  may  be  in  the  form  of  a  specified  amplitude  as  is 
done  in  Section  3  where  the  derivation  of  approximate  pulse  buckling  formulae  is  discussed,  or 
may  be  the  point  of  strain  rate  reversal  detected  by  equation  (25).  Equations  (23,2-1.15  and 
25)  have  been  investigated  via  a  computer  program  for  various  a/h  ratios,  material  parameters 
and  load  amplitudes.  The  loading  function  used  in  the  code  is  derived  from  reference  [9]  whore 
a  parabolic  perturbation  is  assumed.  This  involves  only  cosine,  and  thus  terms,  and  is  of 
the  form: 

6(0.0)  =  t’o  ( 1  -  4:^(cos0  -  cos3^  +  cos  off  -  ...)]  (26) 

20  7TJ  3J  0 

The  maximum  perturbation  is  5  percent  of  v0  and  the  n  =  l  term  is  irrelevent  to  the  asymmetric 
buckled  shape.  The  an  perturbation  terms  are  then  defined  as  on  =  — - — ,  for  n  — 

3,5.7... 


Modei 

a/h 

\'o(in/sec) 

E/,(psi) 

°m  (psi) 

«cr 

t /(sec ) 

<srr  (sec) 

1 

30 

6,000 

1,000,000 

50.000 

17 

0.00243 

0.00194 

2 

30 

4.000 

1,000,000 

50,000 

17 

0.00175 

0.00175 

3 

30 

2,000 

1.000.000 

50,000 

17 

0.00093 

none 

4 

30 

500 

1.000,000 

50,000 

17 

0.00236 

none 

5 

25 

6.000 

1,000.000 

50.000 

14 

0.002 

0.0018 

6 

20 

6,000 

1,000,000 

50.000 

11 

0.00162 

0.00146 

7 

15 

6.000 

1,000,000 

50.000 

9 

0.00121 

0.00121 

8 

10 

6,000 

1.000.000 

50.000 

6 

0.00081 

0.00081 

9 

30 

6.000 

500.000 

50,000 

24 

0.00243 

0.00194 

10 

30 

6,000 

10.000 

50,000 

165 

0.00243 

0.00194 

11 

30 

6,000 

5,000.000 

50.000 

8 

0.00243 

0.00242 

12 

30 

6,000 

1,000.000 

60,000 

19 

0.00207 

0.00165 

13 

30 

6.000 

1.000,000 

70.000 

20 

0.0018 

0.00144 

Table  1:  Results  of  Plastic  Flow  Theory  for  Various  Parameters 

Figures  3  and  4  illustrate  the  growth  of  various  modes  at  several  values  of  time  during  the 
response  for  models  2  and  9  of  Table  1.  It  can  be  seen  that  the  growth  for  several  modes  in 
the  vicinity  of  ncr  is  exponential  with  time,  with  the  ncr  mode  becoming  dominant.  Table  1 
shows  the  critical  mode,  termination  time,  tj ,  and  the  time  at  which  strain  rate  reversal  occurs. 
7 's/?r,  for  several  shell  parameters.  Model  2  (Figure  3).  had  t$rr  =  tj,  so  that  the  solution 
remains  valid  throughout  the  motion.  In  model  9  (Figure  4),  Tj  was  greater  than  ~rrr'.  as  a 
result,  the  solution  was  no  longer  applicable  because  very  rapid  growth  occurs,  as  can  be  seen 
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Figure  3:  Modal  Growth  for  Model  2 


in  Figure  4.  Trends  from  Table  1  indicate  that  the  critical  mode  increases  with  increasing  a/h. 
decreases  with  increasing  tangent  modulus.  Eh.  increases  with  increasing  flow  stress.  erm,  and 
remains  unchanged  with  variation  in  the  magnitude  of  the  initial  velocity  impulse.  The  time 
of  terminal  motion.  Tj  (the  time  at  which  the  shell  velocity  first  reaches  zero),  is  a  function 
of  the  rate  at  which  the  initial  energy  can  be  absorbed  by  the  shell.  As  the  shell  becomes 
thicker  (lower  a/h  ratio).  Tj  decreases  and  as  the  initial  velocity  (energy  input)  increases.  Tt 
increases.  Abrahamson  and  Goodier  [9]  experimentally  investigated  the  behaviour  of  several 
long  cylinders  and  found  good  agreement  with  the  theoretical  mode  predictions  for  lower  a/h 
ratios. 

Reference  [12]  advances  the  theory  of  reference  [9]  to  include  initial  shape  imperfections  in 
the  shell  and  the  effect  of  increasing  t hickness,  decreasing  radius  and  varying  mateTui!  properties 
on  the  response  of  the  cylinder.  The  radius  is  now  defined  as: 


r  =  a  —  rr  0 

where  w0  is  the  uniform  radial  motion.  This  results  in  a  redefinition  of  the  curvature  as: 


1  1  1  d2w 

«  = - +  —(«’-*-  - 

r  a  a 


i2'  ()02  ' 

An  additional  term  in  the  curvature  to  include  the  initial  imperfections  is  given  as: 

1  02u\ 

7*iWi  +  W] 


(21 


(28) 


(29) 
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Figure  -1:  Modal  Growth  for  Model  9 


so  that  equation  (S)  becomes: 


do  _  i  _  2 
OX  p  a 


30 


Using  equations  (27  to  30)  in  equation  (7)  gives  the  equation  of  motion  as: 


£'/,/  dAw  d2w 


•(- 


+ 


r4  '  OO4  Off2 


S.  d2ir  d2u\ . 

+  +  7777 )  +  (,r>  + 


OS2 


oe2 


=  —  w 


02  ir 

Of2 


(31 


where  h  is  now  defined  as  fr,1. 

The  solution  of  equation  (31 )  follows  that  of  equation  (9)  except  that  it  must  be  evaluated 
numerically  if  variations  in  Ett.o.h  and  r  are  to  be  considered.  This  formulation  has  been  used 
in  reference  [12]  to  establish  critical  velocities  for  which  excessive  buckling  will  not  occur. 

An  independent  development  of  the  plastic  flow  theory  is  given  in  Section  2.3  from  Reference 
[11]  which  encompasses  the  elastic  influence  on  the  plastic  flow  behaviour.  The  important 
conclusion  is  that  as  the  a/h  ratio  becomes  higher,  the  elastic  portion  of  the  response  influences 
the  formation  of  the  critical  mode  and  the  ultimate  motion. 
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2.2  Elastic  Pulse  Buckling 

Goodier  and  Mclvor  [10],  and  Lindberg  [ 2 3.-4] ,  produced  a  theory  of  clastic  pulse  bucklina 
for  rings  and  long  shells  subject  to  perturbed,  axisymmetric.  velocity  pulses,  similar  to  the 
perturbation  theory  for  plastic  flow  pulse  buckling.  1  he  mode  formation  occurs  during  elastic 
motion  and  the  theory  considers  only  elastic  behaviour.  Flexural  buckling  modes  are  assumed 
to  form  before  the  hoop  mode  reaches  yield.  Additional  energy  above  that  which  could  be 
absorbed  elastically  in  the  hoop  mode  is  assumed  to  go  into  forming  permanent  plastic  hinsres 
in  the  fiexural  modes.  For  engineering  materials,  the  elastic  buckled  form  can  occur  for  shells  of 
a/h  ratios  greater  than  approximately  2(>0.  Shells  with  an  a/h  ratio  below  this  value  will  have 
their  hoop  mode  enter  the  plastic  flow  region  before  significant  elastic  flexural  motion  occurs, 
although  the  elastic  motion  will  influence  the  resulting  plastic  flow  buckling  mode  [11]. 

The  mechanism  of  producing  buckline  modes  is  a  transfer  of  energy  from  a  fundamental  hoop 
mode  to  flexural  modes.  1  his  is  also  the  physical  mechanism  for  plastic  flow  buckline.  However, 
in  the  plastic  flow,  buckling  theory,  the  hoop  mode  membrane  energy  is  considerably  greater  than 
the  energy  which  is  transferred  to  flexural  modes  and  the  actual  energy  transfer  is  not  considered 
(ie.  the  loss  of  energy  from  the  hoop  mode  is  neglected).  In  elastic  pulse  buckling,  almost  all 
energy  ran  be  transferred  between  the  hoop  and  flexural  modes.  The  differential  equation  of 
motion  for  this  case  is  the  Mathieu  equation  [3j.  which  is  discussed  further  in  Appendix  IF  The 
Ma’hieu  equation  models  transfer  of  energy  from  a  periodic  loading  function  to  a  vibratory- 
buckling  mode  of  one  half  the  frequency  of  the  driving  force.  Parametric  instability  results 
if  sufficient  energy  is  transferred  to  cause  buckling.  Since  we  are  dealing  with  pulse  loading, 
there  is  no  periodic  external  driving  frequency:  however,  the  hoop  mode  oscillates  at  a  natural 
frequency  and  energy  can  be  transferred  to  a  flexural  mode  of  half  of  this  hoop  mode  frequency. 

1  lus  is  termed  ‘autoparametric'  instability  since  the  driving  force  is  self  generated.  The  main 
difference  with  th<>  classical  parametric  instability  is  that  the  autoparametric  instability  case 
does  not  have  a  sustained  energy  input  and  if  the  shell  does  not  buckle  permanently  on  the  first 
phase  of  the  hoop  motion,  it  will  only  vibrate  with  a  continuous  energy  interchange  between 
hoop  and  flexural  modes  until  damping  ceases  the  motion.  The  magnitude  of  the  initial  velocity 
in  relation  to  the  shell  geometry  dictates  whether  permanent  buckling  or  only  vibrations  will 
occur  F.ven  if  buckling  does  not  occur,  stress  values  in  excess  of  hoop  inode  stress  values  are 
generated  by  the  superposition  of  the  amplified  flexural  modes. 

Since  th<'  r>  iipling  between  hoop  and  flexural  modes  is  an  integral  part  of  the  motion,  one 
mode  cannot  be  assumed  to  occur  independently  of  the  other,  as  in  the  plastic  flow  buckling 
theory.  In  plastic  flow  buckling,  the  hoop  mode  thrust  is  assumed  to  remain  constant  until 
motion  stops,  as  it  is  of  considerable  greater  energy  than  the  flexural  modes.  In  the  elastic 
case,  the  transfer  of  energy  needs  to  be  considered  and  coupling  terms  between  the  fundamental 
hoop  mode  and  flexural  modes  have  to  be  maintained  in  the  solution. 
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2.2.1  Derivation  of  the  Equation  of  Motion 

The  equation  of  motion  in  section  2.i  for  plastic  flow  buckling  was  derived  on  equilibrium 
principles.  For  elastic  pulse  buckling,  energy  formulations  and  Lagrange's  equation  are  used  to 
derive  the  equation  of  motion  [Ij. 

The  radial  and  tangential  shell  displacements  with  initial  imperfections,  are  expressed  as: 

TrifU)  =  w(0.t)  +  u-,{0) 

r(0.t)  =  r(0.t)  +  v,(0)  1 32  i 


where  u\  and  e,  are  initial  imperfections  and  ir  and  r  are  the  radial  and  tangential  displacements 
measured  from  the  initial  w,,v,  state. 

The  kinetic  energy  per  unit  length  of  shell  is  given  by: 


Reference  [T  uses  the  dimensionless  quantities: 


_  u  _  u\  _  v  _  c,  _  rt 

a  1  <:  a  '  a  a 


where  c  =  yj' -j-  and  E;  =  r~j.  to  give  equation  (33)  as: 


(33; 


(31 ) 


T  =  -I-iha 


l 


—  )2]cl0. 

'  t it  ‘ 


The  strain  energy  of  the  shell  is  defined  as: 


( 3- Vi 


r  = 


i  30 . 


where  the  hoop  stress  is  given  by  E;(. 

Lite  total  circumferential  strain,  c.  is  the  sum  of  the  middle  surface  hoop  strain.  c ,  and 
the  flexural  strain  from  the  change  in  curvature  a  distance,  r.  from  the  middle  plane: 


(  —  (n.  +  ~A> 

which  when  substituted  into  equation  (30)  gives  the  strain  energy  as: 

1  ,2t 

r  =  -Eihc;  f  [rjj,  -F  cFa'V^Jr/# 

2  Ji) 


(3^1 


where  n2  =  and  K-  the  curvature,  is  defined  by  equation  (2)  and  derived  in  Appendi--  A. 

Equation  (3S)  is  nonlinear  and  terms  up  to  1th  order  are  required  to  maintain  the  required 
coupling  of  hoop  and  flexural  modes  [-1] . 

The  midsurface  hoop  strain  is  defined  by  a  change  in  length  of  a  circumferential  elemenl 
from  add  to  rdd  given  by: 

1  .  Or  „  do  „ 

=  a^\ltdi))  +r  11  m 

which,  when  put  in  dimensionless  form,  and  using  the  relations  w  =  a  -  r  and  v  =  o-  0.  gives, 
after  reduction  [4]: 

dv  do  1  du  ,2 

(,n  ~  ~d6  ~  u  ~  ulje  +  2(dV)  !r 

Including  initial  imperfections  in  the  midsurface  strain  gives:  7rn  =  (,n  -  rm.  wliich  after  sim¬ 
plification  becomes: 


-  (2l 

tm  ~  de 


.do  1  du  >  du  du,  dv  .dv, 

u  ~  ^~do  +  2(tw  *  ~de  ~oe  ~  Ul~de  ~  s~do 


Substituting  equations  (2)  and  (41)  into  equation  (3s)  gives  the  strain  energy  a.-: 

J.  1  r,  f2'r.  do  2  ,#<►’  2  ~.CH-  dudu , 

,  dv  n  dv,  ,  dv  1, 6»u  o  p  o,d2u  0  ,n 

~u'oi  -  -  ?Wrf  +  ‘••'flF*  •*"  l42‘ 

with  shape  imperfections  and  4th  order  terms  maintained.  The  displacements  are  now  expressed 
as  Fourier  series: 

CC 

u{r.$)  =  ci0(r)  +  X>  „(  r )  cos  nd  4  b„(  t  )  sin  nd) 

T1  —  1 

CC 

i'(t.O)  ~  yjcn(T)cos  nd  4  dn(r)  sin  nd]  ( 43  l 

n=  1 

Shell  inextensibility  is  assumed,  which  gives:  cr.  =  —bti/n  and  d„  =  anjn.  Imperfections  are 
taken  in  the  form  [4 1 : 

Of 

—  «,  =  T.  l^n  COS  770  4-  sill  7)0]  (  44  I 


to  give  the  energy  expressions: 


...  .  >  i  ^ ^  It  ■  1  .  l/U,.  > 

I  =  7r/-.//ml( -  Y^(  —  )((— — )‘4(  — 

dr  2  '  iv  dr  dr 


and: 


1  ^ 


V  =  TEiha[al  +  -  ^[(/i2  -  1  )2o2  -  ( n 2  -  2)a0](«?i  -  b'2r 


|  OC‘  oc 

“5  X!  «(n4  ~  3n2  +  4K«n  T  b*)  -  ^(n2  -  2)(*nfln  +  '>„6ri)a0] 

n  =  2 


2S« 


(16) 


To  obtain  the  equation  of  motion,  equations  (45)  and  (46)  are  substituted  into  the  Lagrange 
equation: 

d  or  dr  or 

dt ( di , ]  ~  dsx  +  dx,  ~  0  ('1,] 


to  give: 


and 


1 


ci‘o  +  a0  -  -(»  ~  2)«„  =  0.  for  the  generalized  coordinate,  t,  =  a0 

4 

dn  +  (u-2  -  fnao)an  +  ^n2fng„a*  =  v/12o/ri(!iuao.  for.  i,  =  or, 


(48) 


(49; 


t  2  n2 (n2  —  ]  )2  a2  r  n  2  ( n  2  —  2)  rt  —  T r*  4-  4  ,  »  j  •  •  ,  .  r 

where  u>*  =  — '■  r.:.  +  1 - ,/n  =  — Egn  =  "^(^-2)  ■  anc‘  t'10  nondimensional  imperfection 

is  fn  —  ffn-  The  bu  terms  have  been  omitted  since  solution  will  concern  only  a  few  specific 
amplified  harmonics  and  the  phasing  of  several  harmonic  components  will  not  greatly  affect  the 
solution. 


2.2.2  Impulse  Velocities  for  Vibration  and  No  Buckling 


For  velocities  which  do  not  cause  the  shell  to  buckle  in  the  first  phase  of  the  hoop  mode, 
the  shell  will  vibrate  with  energy  interchanging  between  the  hoop  and  flexural  modes.  At  early 
motion,  most  of  the  energy  will  be  in  the  hoop  mode  and  the  a2  term  can  be  neglected  in 
equation  (T')  to  give  solution: 


I  (i  . 

do  =  • —  Sill  T 
c 


(50) 


for  an  initial  axisymmetric  radial  velocity  of  Vq. 

Substituting  equation  (50)  into  equation  (49)  and  redefining  parameters  in  the  Mathieu 
equation  format  [dj.  gives: 


(in  +  (On  -  /in  sin  r)an  =  VT2 abri  sin  r 
where  the  Mathieu  stability  parameters  are  defined  as  [4j: 


«n  =  <. 


a  v 


( n2  -  n2 


V  +  1  ) 


(51) 


(52) 
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Figure  5:  Mathieu  Stability  Diagram  (from  reference  [10]) 


and: 

W o  n2(n2  -  2)Vq 

c  (n2  +  l)c 

The  format  of  the  Mathieu  equation,  (51),  is  given  in  a  simpler  homogeneous  form  in  refer¬ 
ence  [  1 0  ,  where  initial  imperfections  and  fourth  order  terms  are  omitted  in  the  strain  energy 
expression  since  only  small  velocities  are  considered. 

A  curve  of  Qn  versus  pin,  both  of  which  are  functions  of  the  shell  parameter,  a,  the  initial 
velocity,  Vq  and  the  circumferential  wave  number,  n,  can  be  plotted  on  the  Mathieu  stability 
diagram.  This  is  shown  in  Figure  5.  The  derivation  of  the  stability  regions  (shaded  areas)  on 
the  Mathieu  diagram  is  a  very  difficult  process  and  is  discussed  further  in  Appendix  B  and 
reference  14  . 

It  can  be  seen  in  Figure  5  that  for  this  case  of  small  initial  velocity,  the  curve  passes  mainly 
through  regions  of  stability  with  the  exception  of  the  f2„  =  1/4  abscissa  which  is  unstable  for 
any  nn  (with  no  damping).  Setting  =  1/4  in  equation  (52)  gives  an  approximation  for  the 
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Figure  6:  Response  Showing  Transfer  of  Energy  From  Hoop  to  Flexural  Modes  (from  reference 

W) 


critical  mode  as  [  10] : 


ncr 


1316\/f 


(54) 


assuming  ncr  >>  1.  This  also  gives  the  critical  frequency  of  wn  equal  to  one  half  of  the  hoop 
mode  frequency  necessary  to  have  ‘autoparametric’  instability.  Using  the  relationship  ncr  >>  1, 
the  Mathieu  equation  (51)  can  be  rewritten  in  the  form  [4]: 


rin  -1-  p‘A2(A2  -  sin  r)an  =  p2 A2Q\/l2(5r,  sin  r  (55) 

with  p  —  7^ .  A  =  ~  and  s2  =  ^  from  which  p  can  be  used  as  a  ‘stability’  parameter  in  that  it 
contains  the  necessary  parameters,  a  and  V'0,  to  establish  instability  from  the  Mathieu  diagram. 

The  motion  of  the  shell  is  determined  by  including  the  a2  term  in  the  equation  of  motion, 
(48),  to  allow  for  coupling  between  the  hoop  and  flexural  modes.  Reference  [4]  omits  the 
term  in  equation  (49)  and  attains  a  solution  to  the  shell  motion  by  numerical  integration.  The 
results  are  reproduced  in  Figure  6.  It  can  be  seen  that  energy  is  almost  completely  transferred 
between  the  hoop  and  flexural  modes. 

A  third  curve  in  Figure  6  is  that  of  the  peak  outer  fiber  strain.  This  is  the  structurally 
important  aspect  of  the  autoparametric  behaviour  for  small  initial  velocities  where  permanent 
buckling  does  not  occur.  It  can  be  seen  that  the  peak  strain  occuring  during  the  flexural 
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mode  is  considerably  greater  than  that  during  the  hoop  inode.  An  analysis  not  considering  this 
autoparametric  effect  would  not  detect  this  amplification  and  would  underestimate  the  peak 
stress  under  an  impulsive  load. 

To  estimate  the  stress  amplification  from  flexural  vibrations,  the  peak  hoop  and  flexural 
stresses  can  be  derived  assuming  that  there  is  complete  energy  transfer  uniquely  to  both  modes 
from  the  initial  kinetic  energy.  The  initial  kinetic  energy  can  be  determined  from  equation  (33) 
for  an  initial  velocity.  V'o,  as: 

T  =  irphaVQ  (50) 

The  strain  energy  of  the  shell,  assuming  that  it  is  entirely  in  the  hoop  mode,  is  determined 
from  equation  (46)  as: 

l:H  -  xEihaal  (57) 

The  strain  energy  assuming  that  it  is  entirely  in  the  flexural  mode.  ncr.  is  also  determined  from 
equation  (46)  with  the  additional  assumptions  that  ncr  >>  1  and  that  4th  order  terms  ar<' 
negligible,  as: 

Uj  =  EihanACTa2  a2n  (58' 

Equating  equation  (56)  with  equations  (57)  and  (58)  and  using  equation  (52)  with  =  1/4 
to  establish  the  critical  mode,  gives  an  approximation  of  maximum  amplitudes  of  vibration  as: 


fi0  =  —  and  ancr  =  i/S  — 


(59) 


for  the  hoop  and  flexural  modes,  respectively.  This  means  that  the  flexural  mode  amplitude  is 
2.83  times  the  hoop  mode  amplitude. 

To  compare  stress  amplitudes,  the  maximum  hoop  strain  of  w/a  —  a0  =  V'0 /c  gives  the 
maximum  hoop  stress  as: 

vh  -  Ei—  (60) 

c 

The  flexural  strain  is  defined  from  equation  (37)  for  the  curvature  as: 

h  1  02  w  /'  9  r  i 

<!  =  =  7^ri  Qncr  =  v3on  .r«n  .r  (61) 


which  after  substitution  of  equation  (59)  and  equation  (52)  for  ft,lcr  —  1/4  gives  the  maximum 
stress  in  the  flexural  mode  as: 

<jj  =  v/O  ~Ei  =  VW  (62) 

which  means  that  the  stress  in  the  autoparametrically  excited  flexural  mode  is  2.45  times  that 
of  the  hoop  mode. 
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Goodier  and  Me  Ivor  [10]  numerically  integrated  t  li  e  Mathieu  equation  of  motion,  equa¬ 
tion  (55),  for  several  a/li  values  and  present  graphical  results  which  illustrate  the  energy  trans¬ 
fer  between  hoop  and  flexural  modes.  The  a/h  values  used  are  much  lower  than  those  for  which 
elastic  dynamic  buckling  could  occur. 

The  previous  determination  of  autoparametric  response,  and  of  stress  amplification  factors, 
has  been  for  initial  velocities  such  that  the  a ^  term  could  be  neglected  from  equation  (-19)  for 
the  flexural  terms  of  motion.  Lindherg  and  Florence  [-1]  show  that  for  values  of  the  stability 
parameter,  p  >  1/2,  this  leads  to  unbounded  growth  of  displacement  terms.  The  motion  for 
velocities  resulting  in  p  >  1/2,  can  be  determined  by  numerical  integration  of  the  equations  of 
motion,  (48)  and  (49).  The  effect  of  damping  can  also  be  considered  by  including  3,  the  critical 
damping  parameter,  to  give  the  equations  of  motion  as: 

flj  -f  2 J«0  +  (io  —  -yin2  ~  2)u2  +  2i3<2q  =  u  (Go  ] 

and 

dn  4"  23lVnO-rt  4"  (^'n  fn^o)^n  4”  “H  fnPn^n  —  V^l  2d  fn  t  n  CL  q  (  64  J 

From  results  of  numerical  integration,  it  was  found  that  the  critical  mode  number  was  a  function 
of  p  (ie.  initial  velocity  for  constant  shell  parameters).  An  empirical  lit  to  results  of  the 
numerical  integration  gave  a  formula  for  the  critical  mode  as  [4] : 

«4r  =  ~  +  (0.6)V  (05) 


where,  from  equation  (55),  s2  =  *•. 

From  results  of  the  numerical  integration,  motions  for  two  different  values  of  p  are  shown 
in  Figure  7.  It  can  be  seen  that  energy  transfer  to  the  flexural  mode  is  complete  closer  to  the 
first  inward  hoop  motion  as  p  increases. 

To  derive  an  expression  for  stress  amplification  similar  to  equation  (62)  for  the  small  initial 
velocity  case,  the  total  energy  of  the  system  including  the  4th  order.  terms,  is  used  [4]: 

r^,T^  =  2«o  +  «2  +  2au  +  (7?4°2  -  7,7  (‘0  )ai  +  t~7 77  4  ^  n  =  2p2a2  (66) 

7r  Lina  lb 

In  determining  the  maximum  flexural  stress,  it  is  also  necessary  *o  consider  the  coincident 
energy  in  the  hoop  mode,  arising  from  the  coupling  term,  (n4a2  -  n2ao)a2.  in  equation  (66). 
To  determine  the  maximum  amplitude  of  a o  at  the  point  of  maximum  an.  the  velocities,  oq 
and  dn  are  set  to  zero  and  «2  is  maximized  with  respect  to  (to-  The  upper  bound  to  the  stress 
amplification  is  given  as: 


Of 

—  ~  ^ncr  + 
Oil 


-,  where.  ,4„rr  =  2 


A\4  +  2  -  2A2 
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Figure  7:  Elastic  Motion  for  Various  Stability  Parameters,  p  (from  reference  [4]) 


which  is  shown  to  have  values  in  the  range  of  2.5.  Results  from  numerical  integration  gave 
stress  amplification  values  in  the  order  of  1.65  for  p  >  1.5  indicating  that  the  upper  bound  of 
equation  (67)  may  be  too  high  [4]. 


2.2.3  Pulse  Buckling  from  Large  Initial  Velocities 

For  larger  initial  velocities,  Vo,  in  relation  to  the  shell  parameter,  a,  (ie.  large  p),  the 
shell  will  buckle  during  the  first  compressive  hoop  mode.  For  large  values  of  p,  the  curves  of 
H  versus  p  on  the  Mathieu  diagram  (Figure  8)  do  not  pass  primarily  through  stable  regions 
as  was  the  case  for  smaller  p  values,  as  was  shown  in  Figure  5.  As  was  the  case  in  plastic 
flow  buckling,  many  modes  are  amplified,  with  one  mode  being  dominant.  In  considering  the 
Mathieu  equation,  (55),  it  is  clear  that  the  coefficient  of  an  must  be  negative  for  unbounded 
growth.  Following  the  approach  of  plastic  flow  buckling,  the  largest  n  for  which  unbounded 
growth  may  occur  can  be  estimated  by  taking  sin  r  =  1  in  equation  (55),  to  give: 


A2  <  1, 


V0 


1  2 
nmax  —  S  —  o 
CCt* 


(68) 


Reference  [13]  rewrites  the  Mathieu  equation,  (55),  without  initial  shape  imperfections  as: 


d\  -  p2A2(sin  t  -  A2)a*  =  0 


(69) 
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Figure  8:  Mathieu  Stability  Diagram  for  Large  Initial  Velocities  Causing  Buckling  (from  refer¬ 
ence  [13!) 

with  initial  conditions: 

dn( 0)  =  — ,  an(0)  =  0  (70) 

c 

where  7-  are  harmonic  perturbations  of  the  initial  velocity  field,  Vp.  The  investigation  of  modal 
growth  can  be  accomplished  by  numerical  integration  of  equation  (69)  with  variations  in  the 
parameters  A  and  p.  Figure  9  shows  the  mode  amplifications  with  time  of  an  aluminum  cylinder 
of  a/h=480  and  V0  —  800in/sec  [13].  It  should  be  noted  that  buckling  occurs  for  the  elastic 
case  at  very  high  wave  numbers,  in  the  range  of  n  =45  to  65. 

For  complete  analysis  of  the  motion  of  the  elastic  shells,  equations  (48)  and  (49)  need  to 
be  numerically  integrated  with  the  initial  conditions  of  equations  (70).  An  empirical  fit  to 
results  obtained  in  this  manner  gives  a  relation  for  the  amplitude  of  the  most  amplified  mode 
in  relation  to  the  initial  imperfection  as  a  function  of  the  parameter,  p,  as  [4j: 

Y-\max  =  1.2e”  (71) 

On 

The  mode  of  maximum  amplification  can  be  estimated  by  solving  equation  (69)  for  a  con¬ 
stant  hoop  thrust  case  (ie.  letting  sin  r  =  1).  This  is  the  same  assumption  as  is  made  in 
plastic  flow  buckling,  where  the  time  dependent  coupling  between  hoop  and  flexural  modes  is 
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Figure  9:  Modal  Amplification  for  Elastic  Dynamic  Buckling  of  Thin  Walled  Cylinder  (from 
reference  [13]) 


neglected.  This  gives  [13]: 


QA(f) 


1  oAn  sinli  Tnr 


where  Tn  =  pX\/\  —  A2 


(72) 


from  which  the  maximum  amplified  mode  can  be  determined  by  minimizing  equation  (72)  with 
respect  to  A  to  give: 


A  cr  — 


v/2 


(73) 


To  establish  whether  elastic,  dynamic  buckling  will  occur  for  a  given  velocity  pulse  on  a 
given  shell,  it  is  necessary  to  determine  if  either  hoop  or  flexural  yielding  occur  and  which  will 
occur  first.  Elastic,  dynamic  buckling  is  signified  if  flexural  yield  is  reached  first,  otherwise  the 
shell  will  enter  the  plastic  flow  hoop  mode  where  dynamic  buckling  will  occur  as  plastic  flow 
buckling  influenced  by  the  elastic  mode.  L  sing  equation  (61 )  for  flexural  strain  and  equation!  71 ) 
for  maximum  amplitude,  the  flexural  strain  can  be  expressed  as: 

c/  =  \/3A2p(1.2Vp)J  (74) 


Setting  expression  (74)  equal  to  the  yield  strain  cy  gives  a  relation  between  p  and  h/a.  A  similar 
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Figure  10:  Hoop  and  Flexural  Yield  and  Dynamic  Buckling  Limits  for  Elastic  Response  to 
Pulse  Loads  (from  reference  [4]) 


expression  for  the  hoop  strain  can  be  attained  by  setting: 

Vo  Ph 

—  =  ^J=~  =  £y  (  <  5) 

C  v/l2a 

Using  the  value  of  Acr  =  0.6,  from  results  of  numerical  integration  of  equation  (69),  Lindberg 
and  Florence  [4!  present  a  curve,  Figure  10,  for  different  imperfection  values,  6„,  giving  hoop 
and  flexural  yield  limits  as  functions  of  p  and  a.  For  practical  values  of  initial  imperfection 
of  1  percent  of  the  shell  thickness  (<5  =  .01)  and  ev  =  .004,  it  can  be  seen  that  purely  elastic, 
dynamic  buckling  will  not  occur  for  shells  below  a/h=260.  The  stability  parameter,  p,  must  also 
be  greater  than  3.66  for  buckling  to  occur,  otherwise  only  autoparametric  vibrations  will  occur. 
These  parameters  represent  a  very  thin  shell,  thinner  than  would  be  used  in  most  submarine 
applications. 

The  topic  of  autoparametric,  dynamic,  elastic  buckling  has  been  addressed  subsequent  to 
reference  [10]  by  other  authors  Hubka  [1],  derives  an  analytical  solution  to  the  coupled  non- 


linear  differential  equations  of  motion  for  response  to  small  initial  velocities.  The  solution 
agrees  with  the  numerical  results  of  reference  [10]  showing  a  slowly  varying  hoop  mode  function 
modulating  a  rapidly  varying  flexural  mode. 

The  Mathieu  type  instability  describes  response  for  periodic  type  loadings  and  many  ref¬ 
erences  examine  this  [15,16,17].  This  is  not  pursued  further  here  as  response  to  single  pulse 
loading  is  of  primary  interest. 
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2.3  Elastic  and  Plastic  Pulse  Buckling 

This  section  presents  an  analytical  investigation  [11]  which  includes  both  the  plastic  and 
elastic  regimes  in  a  single  theory.  The  theory  reduces  to  the  plastic  flow  or  elastic  theories 
discussed  in  the  previous  two  sect  ions  for  limiting  cases.  The  necessity  of  combining  both  elastic 
and  plastic  behaviour  for  a  wide  range  of  a/h  ratios  is  demonstrated.  Although  equations  of 
motion  are  produced,  a  general  solution  to  the  elastic-plastic  case  is  not  formulated. 

The  problem  considered  here  is  again  for  a  ring  or  infinite  shell  subject  to  a  perturbed, 
nearly  uniform,  axisyrnmetric  velocity  pulse.  The  linearized  equation  of  motion,  derived  from 
equilibrium  of  a  shell  segment  (sec  equation  7)  is.  with  the  additional  assumption  that  u  is 
much  less  than  1: 

ii  a  /  E 


U'V  +  (1  +  + 


+ 


=  0 


(76) 


a232  a232  '  o2  32 

where  32  —  and  is  the  spatial  derivative  with  respect  to  0.  The  parameter  J.  in  this  form 
of  the  equation  of  motion,  can  be  used  for  the  elastic  as  well  as  the  plastic  flow  case.  For  purely 
linear,  elastic  response.  3  =  1  and  o/E  -  c.  Reference  [11]  considers  only  material  of  bilinear 
stress-strain  curves  represented  by: 


cr / E  =  c,  32  =  1,  for  0  <  c  <  (y 

a/ E  =  (1  -  32) ty  +  J-t.  32  =  for  cy  <  c 

Equation  (76)  is  then  divided  into  two  regimes,  one  for  elastic  strain: 

u"  +  ( 1  -r  —  )u"  +  A;  +  — r  =  0.  for  0  <  c  <  cv 
q*  q*-  a2 

and  one  for  plastic  strain: 


(77) 


I'M 


(1 


1  -  T2  c, 


32  o  - 


+ 


a232 


+ 


1  ~  32  c, 

32 


^  =  0.  for  cv  <  c 


(79) 


The  initial  conditions,  as  in  the  plastic  and  elastic  cas,5s.  incorporate  imperfections  in  the  initial 
displacement  and  velocity  fields,  defined  by  t he  series: 

l'o 


tivi  8)  —  — -  Y(atA)  sin  nO  —  hT|0  cos  n9) 


Vo 


n=i 

oe 


u0(6)  =  —  [1  +  y>n,  sin  nff  -f  bn,  cos  n^)j 


(80) 


n  =  1 


where  on,  are  of  perturbation  order  (<<1). 
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TV  solution  follows  t hat  of  the  previous  rases,  with  a  particular,  unperturbed  solution  and 
a  perturbed,  complimentary  solution  yielding  buckling  terms. 

Solving  equations  (78)  and  (79)  for  perfectly  radial  motion  (no  dependence  on  6).  yields  the 
particular  solutions: 

_  1  o  .  ,  1  o  .  ,  .  . 

u=  — si  nr.  for  —  =  u  <  cv  (81) 

c  c 

and: 

1  -  J2  .  <y  /l  - 


u  - - cv  +  “7 


J2  'y 


3  V  3 2 


—  +  ( - )*  sin(  3(  r  -  r„)  +  ^i.  for 


where  rf)  =  arcsin  cot  ^  =  3y(-^-)2  —  1  and  the  constants  of  integration  were  derived 
from  the  initial  conditions: 

>  1  o  •  ,  —  ,  1  o  ,  „ 

u(tq)  =  — sin  r0  =  c. .  and  ,u(r0t=  —  cos  Tp  (  8,5 1 

c  c 

If  the  initial  velocity  is  great  enough  to  cause  \  o/c  >  <v.  then  the  radial  hoop  motion  ceases  at 
time: 

jt/2  -  r 

r  j  =  r o-l - — - .  (s-l) 

The  complimentary  solutions  are  obtained  by  solving  equations  (78)  and  (79!  subject  to 
initial  conditions  (SO).  The  flexural  motion  of  the  shell  is  assumed  to  be  ine.xtensional.  allowing 
(  =  u.  and  the  complete  solution  to  be  taken  as  u[9.t )  =  t7(r)  -f  uc(9.r).  Equations  ^78)  and 
(79).  reduce  to  the  variational  equations  of  the  perturbation  of  17: 

u'3  4-  ( 1  +  — rr  )u  4-  —  =  0  for  0  <  u  <  cy  (8.r>! 

Q  2  a  - 


1  “  (V  ,  ,  «  , 

o-”“ 


c—  =  0  for  <  v  <  77 


for  which  solutions  of  the  form: 


i,.( 9 .  r )  =  —  Y"  [au'  r  )  sin  n9  +  brJ  r )  cos  n9] 
c 


are  assumed. 

Substitution  of  equation  (87)  into  equations  (85)  and  (86)  produces  the  homogenc  ms  form 
of  the  Mathieu  differential  equal  ion .( 51 !.  of  the  e!a>lic  theory  (section  2.9).  if  n<<l: 


f/r  4-  <>(n4  -  n2  — ,^i n  r)r/i;  =  0,  for  0  <  u  <  cv 
a  - 


and  of  the  plastic  regime: 


gn  +  Q2J2{n4  +  1-_L(  _^_)2  sin[J(r  - 


To)  T  Vj}<;„  =  o.  for  <y  <  u  (89i 


the  latter  of  which  reduces  to  : 


a2Eh 


—  \gn  =  0.  for  i2  <<  ] 


(90) 


Equation  (90)  is  similar  to  equation  (ID)  of  the  plastic  How  theory  (section  2.1  h  Equation  •  sf> ) 
reduces  to  equation  (8-8)  for  elastic  motion  where  ,3=1. 

Both  equation  (88)  for  the  elastic  case  and  equation  (89)  for  the  plastic  flow  case,  are 
of  the  Mathieu  instability  type.  Values  of  the  parameters,  n.o2.V0.c  and  J  can  be  found 
which  fall  within  the  unstable  regions  of  the  Mathieu  diagram.  If  the  hoop  mode  response 
remains  completely  elastic  (ie.  r0  <  c<v),  then  the  mode  of  buckling  will  be  established  solely 
by  equation  (88)  as  was  given  in  section  2.2  from  Lindberg  and  Florence  [4],  If  \'u  >  c( .  . 
then  the  hoop  mode  enters  the  plastic  material  regime  and  the  buckling  mode  shape  will  be 
some  combination  of  the  two  predominant  harmonics  determined  from  the  elastic  and  plastic 
solutions.  If  Vq  >>  c(y,  then  the  response  will  be  governed  by  the  predominant  mode  of  the 
plastic  flow  regime  as  was  the  case  in  section  2.1. 

To  obtain  the  predominant  harmonics  of  equations  ( S S )  and  (S9).  these  Mathieu  equations 
must  be  solved.  This  is  of  extreme  analytical  complexity  [14].  Stuiver  [11],  produces  an  ap¬ 
proximate  solution  resulting  in: 


Aei  = 


n .  .• 

s 


for  the  elastic  regime,  where  A .p  and  s  are  defined  in  equation  (55)  and: 


(91) 


’  ^  fX i .  1  /o  , 

;  =  0.85vGr(-f  )'A  1  +  (3 

It  t 


l'o/e 


2’  1  /■! 


i  92  | 


for  the  plastic  regime.  lor  p  >  2.  Ae;  is  equal  to  0.595  which  is  in  good  agreement  with  the 
value  of  Af/  =  0.6  obtained  through  numerical  integration  [4],  Equation  (92).  for  the  plastic 
flow  wave  number.  n;j;.  gives  comparable  values  to  equation  (24)  of  the  plastic  flow  case  derived 
by  Abrahamson  and  Goodier  [91. 

Comparison  of  resultant  buckled  wave  numbers  from  tests  on  cylinders  with  9  <  ajh  <  36 
show  that  the  experimental  values  lie  between  the  nfi  and  npj  values  of  equations  (91)  and  (92) 
with  no  particularly  good  agreement  with  either  value.  Stuiver  postulates  that  this  is  a  result 
of  both  the  elastic  and  plastic  modes  influencing  the  response. 
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A  complete  solution  to  the  nonlinear  differential  equations  describing  the  radial  motion 
would  be  very  complex  and  would  be  better  left  to  approximate  numerical  methods,  such  as 
finite  difference  or  finite  element. 

One  such  study,  using  a  dynamic  elastic-plastic,  geometrically  nonlinear  finite  difference 
code  for  rings,  is  presented  by  Wesenberg  [IS].  Three  cylindrical  shells  of  a  h=100.  200  and 
300  were  subjected  to  axisymmetric  impulse  loads,  and  response  was  measured  by  high  speed, 
photography.  Initial  imperfections  were  measured  for  the  predominant  resultant  mode  and 
incorporated  in  the  finite  difference  solution.  The  radial  velocity  perturbations  were  measured 
from  the  photographs  and  used  as  the  velocity  perturbation  parameters  in  the  finite  difference 
solution.  The  numerical  solution  allowed  the  shells  to  buckle,  deform  plastically  and  damp)  out 
to  a  final  deformed  shape.  The  permanent  peak  deformations  for  a  range  of  impulse  loads  were 
in  good  agreement  with  the  finite  difference  predictions.  The  wave  number  of  predominant 
response  was  shown  to  increase  linearly  with  the  magnitude  of  the  applied  impulse.  All  three 
shells  buckled  in  wave  numbers  ranging  from  30  <  n  <  35;  however,  the  solution  does  require 
imperfection  and  perturbation  input  from  the  buckled  shell  which  one  would  normally  not  know 
apriori. 

Using  equations  (91)  and  (92)  for  the  experimental  aluminum  shell  described  in  reference 
11.  with  parameters:  h=0.02  in,  a=2.0  in.  E=10, 000.000  psi,  £/,  =  130.000  psi.  <7(0.291 ) =44 .000 
psi.  f>  =  0.27  X  10~3  lb-scr2/tn4,  c=200.000  in/sec  and  Uo=805  in/sec.  gives  nei  =  13  and 
npi  =  121.  This  shell  buckled  experimentally  at  n  =  30  [18].  This  indicates  that  the  actual 
response  is  somewhere  in  between  the  cases  of  purely  elastic  and  purely  plastic  response  for  a 
shell  of  a/h=100.  The  plastic  theory  greatly  overestimates  the  mode  of  fundamental  response. 
Since  both  nf>  and  np;  are  linear  functions  of  a/h,  for  the  same  velocity  impulse,  the  shell  with 
a/h  =  200  will  have  ?<,/  =  20  and  nPi  =  242  as  predictions  and  the  a/h  =  300  shell  will  have 
n,j  --  39  and  n pi  =  363  a-  predictions.  Tor  these  thinner  shells,  which  buckled  experiment a’ly 
at  n  =  32  and  n  =  33.  respectively,  the  elastic  prediction  is  in  reasonable  agreement.  This  could 
be  expected  as  these  shells  are  in  the  range  of  the  elastic  limit  of  a/h>260  reported  by  Lindberg 
and  Florence  [4], 


2.4  Pulse  Buckling  for  Finite  Length  Cylinders 

The  theories  presented  in  t he  previous  sections  of  this  chapter  have  dealt  with  rings  or 
long  shells  for  which  plane  strain  conditions  have  been  assumed.  This  section  discusses  theory 
for  the  dynamic  buckling  of  axisymmetrically  loaded  shells  of  finite  length  where  biaxial  stress 
states  and  end  conditions  have  been  included.  First,  the  extension  of  the  plastic  flow  buckling 
case  to  include  a  biaxial  stress  state  will  he  discussed.  Then,  the  elastic  dynamic  buckling  case 
is  considered  for  finite  length  cylinders  with  simply  supported  end  conditions.  Elastic  dynamic 
buckling,  particularly  if  the  loading  function  is  of  longer  duration  than  a  put’.'  impulse,  may 
occur  at  circumferential  wave  numbers  of  small  enough  value  that  end  conditions  affect  the 
response  of  the  entire  shell  [19]. 


2.4.1  Plastic  Flow  Buckling  of  Finite  Length  Cylinders 


As  was  the  case  in  section  2.1  for  rings  or  long  cylinders,  the  theory  derived  for  finite 
length  cylinders  is  for  the  specific  case  of  a  perturbed,  axisymmetric,  velocity  pulse  on  an 
imperfect  cylinder.  The  theory  is  most  applicable  to  thick  shells  [a/h  <  -40).  where  the  hoop 
membrane  strain  will  undergo  significant  plastic  deformation  before  buckling  occurs  and  elastic 
behaviour  can  be  neglected.  The  assumption  of  no  occurrence  of  strain  rate  reversal  is  also 
made  and  a  linear  strain  hardening  modulus  is  used  as  the  constitutive  relationship.  Florence 
and  Vaughan  [20].  first  developed  the  theory  for  short  shells  with  a  plane  stress  assumption, 
and  then  extended  it.  to  consider  shells  of  variable  length  to  the  limit  of  the  infinite  length, 
plane  strain  case  [21].  The  solutions  are  for  the  unbounded  amplification  with  time,  of  the 
predominant  circumferential  buckling  mode. 

Due  to  the  biaxial  stress  state,  an  additional  resisting  moment  to  the  strain-hardening 
moment  derived  in  section  2.1  exists.  Material  incompressibility  is  assumed  giving  the  condition 
(in  terms  of  strain  rates): 

(x  +  (f.  +  =  0  (9d) 

where  c-  =  The  generalized  strain  rate  and  stress  are  defined  by  [29.22’: 


and: 


/2=  |(<J +  <?  +  <?) 


+  a'x  ) 


(9-n 


(95) 


where  a[  are  the  deviatoric  stress  components.  With  the  assumption  of  n.  =  0  sires.-,  through 
the  shell  thickness),  and  using  the  Levv-Mises  flow  law  [22].  i,/a[  =  A.  where  A  =  3c/2a.  from 
equations  (91)  and  (95).  the  stress  components  become: 


2a  2a  . 

aT  =  — r(2 cr  +  co).  and.  ap  =  —  (Cj-  +  2<p) 
dr  dc 


(90) 


dl 


and  for  a  linear  strain  hardening  assumption,  a  =  which  allows: 


a  -  ay  +  Eh  t 


(07) 


where  ay  is  the  yield  stress. 

For  unperturbed,  radial,  hoop  motion,  the  circumferential  strain  rate  is: 

(0  =  -(1  -  2/fl)(u’o/fl)  (9V) 

where  .  is  measured  from  the  shell  midsurface  and  ti'o  is  the  radial  hoop  motion. 

With  the  assumption  of  plane  stress.  aT  —  az  —  0.  and  (T  =  c:  =  —  tf/2  from  equation  (9(i). 
at  the  midsurface,  equation  (98)  gives  cx  =  c,  =  ii0/'2a.  The  strain  rates  at  any  point  r.  in  the 
shell  are  defined  as: 


(T  =  ti-Q /2a.  (e  =  -(1  -  ~/a )( u'o/a).  and.  e.  =  -(1/2  -  z/a)(i.i:0  a)  (99) 

which,  when  substituted  into  equation  (91),  with  (r/a)  terms  of  order  two  or  more  neglected, 
gives  the  generalized  strain  rate: 


f  =  -(1  -  z/a)(w0/a)  (100) 

and  the  stress  relations: 

<jr  =  ( 2z/3a  )rr  and  <Je  =  -( 1  -  :/'2a  )a  ( 101 ) 

where  the  generalized  stress  is  now  defined  as: 

rr  =  ay  +  £/,(  i  -  c/fl)(u-0/ai  (102) 

Reference  ; 20]  defines  an  as  the  current  inidsurface  stre:.-: 

<7°  =  ery  4-  Eh(  U'o/fl  )  (103) 

allowing  the  stress  components  to  be  defined  as: 

aT  -  (2c/3a)a0  and  o0  =  -<7°  +  {z/'io)(av  +  4£/ltr0/«)  (104) 

again  with  :/a  terms  of  order  two  or  greater  neglected. 

The  moments  and  membrane  forces  for  a  shell  segment,  determined  using  equation  (104) 

are: 

h/2 

Mt=  I  <r-;dz  =  {liJ/lSa)o°  (103) 

J-h/2 
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where  w ,  is  defined  in  reference  [20]  as  the  initial  deformed  shape,  w,  =  u>(#,0).  The  generalized 
stress  is  then  given  as: 


o  ,  Eh,  s  Ehz '  ,  ,  ,  ,  d2(w  -  tij,-), 

a  =  o  +  —(w  -  it-,)  -  -^-[u>0  +  (w-  wt)  +  - £02 - i 


(115) 


Equations  (96),(  112)  and  (115)  and  a  binomial  expansion  of  equation  (113),  give  stress  compo¬ 
nents: 


ax 


eg 


2 z,  0/,  1  d2u>  Eh.,  .. 

—  <r°(l  +  — -r^r)  +  — (tv  -  u\) 
3a  w0  86 2  a 


_[„»  +  -  „,)]  +  £[*°(1  +  ~gijT)  +  T (3“’°  +  4<“’  -  Wi)  +  3  W 


(116) 

<92(u’  -  it-,) 


with  resultant  moments  and  membrane  forces: 


(117) 


+  j-0w] + ^<3u’° + 4(“’  ~  “’■> + 3?(^'",l)>1 


(118) 


Nx  =  0  and  jV$  =  -/i[<r°  +  — (tn  —  «?,•)]  (119) 

a 

The  second  term  of  Mg  in  equation  (118)  is  the  resisting  moment  from  strain  hardening,  Eh, 
which  governed  the  theory  of  section  2.1  for  infinite  shells.  The  first  term  of  Mg,  ^a°(  1  + 
‘s  an  additional  resisting  moment  arising  from  the  different  locations  of  points,  z, 
through  the  shell  thickness  on  the  yield  ellipse  (Figure  11),  which  is  defined  by  the  equation: 


<72  -  azug  +  a g  =  a2 


(120) 


The  strain  rate  vectors  at  the  outer,  mid  and  inner  surfaces  are  of  different  magnitudes  on 
the  yield  surface  giving  rise  to  a  strain  differential  through  the  shell  thickness  and  thus  a 
resisting  moment.  This  moment  has  been  termed  the  ‘directional  moment’  in  reference  [20] 
and  is  dependent  on  a  biaxial  stress  state.  The  second  term  of  Mg  is  referred  to  as  the  ‘strain 
hardening’  moment  and  depends  on  the  existence  of  Eh- 

The  equation  of  motion  is  found  by  substituting  the  curvature  and  force-moment  expressions, 
equations  (111)  and  (117)  to  (119)  into  equation  (7),  giving: 


1  d2Mg 

a 2  d02 


-M  + 


1  d2w .  d2(wo  +  u-) 

^Wl  +  ph  &t2 


=  0 


(121) 
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,Z  --  -  h/2 
>Z  •  0 
tZ  -  h/2 


Figure  11:  Strain  Vectors  on  Yield  Ellipse  (from  reference  [20]) 

Simplifying  assumptions  are  introduced  by  Lindberg  and  Vaughan  [20];  the  hoop  thrust  is  made 
constant  with  time:  Ng  =  —o oh,  the  (w0  +  w)/a2  term  in  the  curvature  is  neglected,  and  the 
4 (w  -  W{ )  term  in  Mg  (equation  (118))  is  neglected.  This  yields  equation  (121)  as: 


ft3  f  0/-t  ,  1  d2w.  Ehd2(w-w,)  <7q  d2w  d2w  _ 

36a['  (1  +  w0d62>  +  3a  362  ]+a2d02+pdt2  ° 


(122) 


where  the  four  terms  represent,  respectively,  the  directional  moment,  the  strain  hardening  mo¬ 
ment,  the  hoop  thrust  and  the  inertia.  With  the  introduction  of  the  dimensionless  parameters: 
ti  =  w/a,Ui  =  w,/a,u0  =  wo/a,T  =  Vot/2a,rj  =  pV^  /2<ry,a2  =  h2/\2a2  and  3  =  Eh/vy 
equation  (122)  becomes: 


Q2  d2  .  1  d2u  p2(u  -  u,)1  T/^U  _ 

c  do2  [uQ  ee2  p  de 2  J  ee2  2  dr 2 


(123) 


which  is  similar  to  equation  (41)  of  the  plane  strain  plastic  flow  buckling  case  except  that  the 
‘directional  moment’  term  has  been  added  and  u  has  been  neglected  with  respect  to  1. 

The  equation  of  motion  for  the  unperturbed  motion,  equation  (108),  is  given,  from  the 
simplifying  assumptions,  as  ay  =  — aptiio,  which,  when  put  in  nondimensional  form  gives:  uq  = 
— 2/t/.  Initial  conditions,^  =  0  and  wq  =  Vo  yield  uq  =  0  and  ito  =  2.  Integration  of  tio  with 
the  initial  conditions  yields  the  dimensionless  hoop  velocity: 


u0(r)  =  2(1  -  r/r/) 


(124) 
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and  the  dimensionless  hoop  displacement: 


«o(r)  =  2 (r  -  t2/2tj). 


(125) 


The  time  at  which  motion  ceases  (uo  =  0),  is  t/  and  the  midsurface  hoop  strain  at  this  time 
is  wo/a  =  'io(Tf)  =  Tf.  Reference  f20l  introduces  a  further  parameter,  £  =  1  -  t/tj ,  s-'h  tKat 
the  period  of  motion  is  now  defined  in  the  interval  0  <  f  <  1.  The  equation  of  motion  is  now 
expressed  as: 

W + 2r‘w  -  £k  -  2^l?<“  -  “•»!  -  »■  <126> 

The  perturbations  of  the  fundamental  motion  are  expressed  in  sine  series  as: 


=  ^u„(£)sinn0 


(127) 


with  initial  shape: 


and  initial  velocity  perturbations: 


«.(<u)  =  E  an  sin  nd 


V'  =  f/0(l  +  £>„  sin  nO) 


(128) 


(129) 


yielding  the  Bessel  type  differential  equation: 


subject  to  initial  conditions: 


..  Qn  .  n2  c 

Un  —  —^~un  —  Rnun  —  Jn0.n 


u„(  1)  =  an  and  ti„(l)  =  -2 rfbn 


where: 


Qn  =  — - — ,  R2  =  2tjti2(1  -  a2(3n2),  and  sn  =  2 rfa2f3n4 


(130) 


(131) 


(132) 


The  solution  of  equation  (130)  takes  the  form  of  Bessel  functions  and  is  [20],  for  R2  >  0: 

«n(0  =  An(Z)an  +  Bn(£)bn 


(133) 


where: 


^n(0  =  2TfC[Tv(RnOKv-l(Rn)  +  Kv(RnOIv-l(Rn))(Rn  +  Sn/Rn )  -  Sn/R2n 
Bn(0  =  2Tfe[-Iv(RnZ)Kv(Rn)  +  Kv(Rnt)Iv(Rn))  (134) 
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Figure  12:  Displacement  and  Velocity  Amplifications  as  a  Function  of  Harmonic  Number  for  A 
Short  Cylinder  (from  reference  [20]) 

where  /„  and  Kv  are  modified  Bessel  functions  of  the  first  and  second  kind  of  order  v  =  (— fr ). 
The  restriction  f?2  >  0  defines  the  range  of  harmonics,  n,  which  will  grow  unbounded,  as 
n2  <  l/a2/3.  For  i?2  <  0,  the  solution  is  in  terms  of  Jv  and  Yv  unmodified  Bessel  functions 
which  do  not  become  unstable.  This  is  similar  to  the  plane  strain  case  where  hyperbolic 
functions  represent  the  unbounded  range  of  n,  and  circular  functions,  the  bounded  n.  The 
An(£)  and  Bn(£)  functions  represent  amplifications  of  the  shape  and  velocity  imperfections, 
respectively.  Equation  (130)  may  also  be  evaluated  numerically.  The  amplification  functions 
A„(£)  and  Bn(£)  resulting  from  numerical  integration,  are  shown  in  Figure  12.  These  are  similar 
in  form  to  Figures  3  and  4.  The  exact  values  of  the  amplified  parameters  are  dependent  on  the 
initial  shape  and  velocity  imperfections.  Although  the  displacement  amplification  functions  are 
greater  than  the  velocity  amplifications,  the  two  terms  may  be  comparable  for  realistic  initial 
imperfections.  An  approximate  formula  for  the  critical  mode  is  given  in  reference  [20]: 


Reference  [21]  carries  the  theory  for  short  shells  to  shells  of  variable  length  by  introducing 
the  parameter,  k,  into  the  strain  rate  relation,  such  that: 

ex  =  -keg  for  0  <  k  <  1/2  (136) 
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Figure  13:  Yield  Ellipse  for  Variable  Length  Shells  (from  reference  [21]) 


where  k=l/2  for  the  plane  stress,  short  shell  case  just  described,  and  k=0  for  the  infinite  shell, 
plane  strain  case.  Figure  13  shows  the  yield  ellipse  of  the  midsurface  strains  for  varying  shell 
length.  The  shaded  region  is  entirely  in  a  state  of  biaxial  compression  with  the  short  shell  case 
(k=l/2),  corresponding  to  ax  —  0.  The  stress  state  is  expressed  in  terms  of  k  as: 


Q j. 

OQ 


1  -  2k 

2  -  k 


(137) 


The  value  of  k  varies  along  the  length  of  the  shell,  being  near  zero  at  the  middle  of  a  long  shell 
and  near  one  half  at  the  ends  of  the  shell.  An  approximate  formulae  for  k  as  a  function  of 
length  is  given  [4],  as: 

, ,  „  cosh(2flf;) 

(138) 


,  ,  .  cosh(2flfj) 

k(x)  =  for  ‘  1,2  s  1  s  1/2 


where  L  is  the  cylinder  length,  D  is  the  diameter  and: 


ft2  = 


(2  -  k)Kx 

[K\  —  (1  —  2A:)2] 


(139) 


The  theory  for  the  variable  length  shells  follows  that  of  the  short  shells,  but  includes  addi¬ 
tional  narameters  as  a  result  of: 

ix  =  kwo/a  (140) 
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which  are  defined  as: 


K\  =  2(1  -  k  +  k2),  K2  =  \/3A'i/2,  and,  K3  =  (2  -k)/Kx  (HI) 

Only  the  key  equations,  with  the  K,  parameters,  will  be  repeated  here  .  The  solution  of  the 
unperturbed  radial  motion  is  delined  by  equation  (110)  which  now  becomes: 

—  =  sin(Ac/,</ a)  -  -—^[1  -  cos(A cht/a)]  (142) 

a  ACh  2h2k'h 


where  A2  =  |(2  —  k).  This  solutionis  reduced  to  that  of  equation  ( 125)  for  a  constant  midsurface 
stress,  om,  giving: 

T/1a 

(143) 


y  =  r( 2  -  T/Tf) 


with  tj  now  defined  as: 


TJ  = 


PvoK- 


2am(2  —  k)  ’ 


(114) 


crm  is  taken  as  the  average  stress  in  the  post  yield  region  as  was  used  in  section  2.1.  This  is 
different  than  the  short  shell  theory  which  uses  <ry,  which  is  the  yield  stress  (ie.<rm  >  ay). 

The  resultant  moments  and  forces  for  the  perturbed  motion,  equations  (117  to  119),  become: 


Mx  = 


h3  3kamK3 
12a[(2-k)K2{  + 


1  d2w 
ii)Q  dd2 


)-2(2k-l) 


K3Eh 
3  a 


Q 2 

-  Vli))] 


(145) 


12a  (2  —  k)Kj'  Wade2 


)  +  2(2  - 


«■•<))]  (146) 


Nx  =  ^(2k  -  l)  B.nd  Ne  =  ^(-2  +  k)  (147) 

A  2  A  2 

where,  again,  Mg  consists  of  a  ‘directional  moment’  resulting  from  a  strain  differential  on  the 
yield  surface  through  the  shell  thickness  and  a  ‘hardening  moment’  resulting  from  an  Eh  of 
positive  slope  creating  a  strain  differential  through  the  shell  thickness. 

The  above  force,  moment  and  curvature  relations  are  substituted  into  the  equation  of  motion 
(equation  121),  resulting  in  the  nondimensional  equation  of  motion: 


d2u  n  d2u  2  d2  .  3k2I\3  cPu  2(2  -  k)r j (3  d2 u  d2U{ 

0E2  +  T}W2  ~  Q  de2\2  -  k)2£.  &td62  ~K2 


(148) 


where  the  nondimensional  quantities  are  defined  in  equation  (123)  except  that  am  replaces 
o y  and  Tj  is  defined  by  equation  (144).  Using  the  same  displacement  perturbations  (equa¬ 
tion  127),  initial  distortion  (equation  128),  velocity  imperfections  (equation  129),  and  initial 


39 


conditions  (equation  131),  the  same  Bessel  type  differential  equation  as  equation  (130)  results. 
The  coefficients  are  now  redefined  in  terms  of  k  as: 


Qn 


3a2nAk2K$ 

(2  -  k? 


Rl  =  2r/n2[l  -  (2  -  k)^p-],  and,  5n  =  2 rfn2  -  R\ 

A  2 


(149) 


with  solution  defined  by  equations  (133)  and  (134).  The  amplification  functions  -4„(£)  and 
Bn(£)  behave  as  in  Figure  12. 

Both  references  [21]  and  [20]  develop  simplified  expressions  for  the  amplification  functions 
by  approximating  their  final  values  at  £  — ►  0.  In  considering  the  modified  Bessel  functions  at 
their  limits,  it  can  be  shown  that: 


Lt^oCh(RnO  =  0  and  Lt^0CKv(Rni)  = 


2v~lr(v) 

Ra 


where  r(t>)  is  the  Gamma  function,  and  for  Rn  >>  v,  it  can  be  shown  that: 


h(Rn)  ^ 


\J2  nRn 

which  result  in  approximations  for  An  and  Bn  at  £  =  0  of: 

.  r(t>)e*n77n2  ,  T(v)eRnTt 

”(0;  “  4 ( 7T ) 1  /2( Rn/2 ) "+3/ 2  an  n(0)  -  2(T)l/2(i?7i/2)v+l/2- 


(150) 


(151) 


(152) 


For  k  approaching  a  value  of  zero  for  the  long  shell,  plant  -train  case.  Qn  =  0,  i?2  = 
2r/n2(l  -  2-^-),  and  5„  =  4  and  the  solution  of  equation  (130)  now  is  in  the  form  of 

hyperbolic  amplification  functions: 


An(0  =  (2tjti2 / A2)  cosh  Rn(  1  -  ()  -  Sn/R2n  (153) 

and: 

Bn(0  =  {2rf  /  Rn)  smh  Rn(l  -  f)  (154) 

which  for  terminal  motion  (£  =  0),  are  approximated  by: 

>ln(0)  ~  TjeRnn2/R2  and  Z?n(0)  TjeRn /Rn  (155) 


Differentiating  either  of  equations  (155),  with  respect  to  n  to  obtain  its  maximum  value  yields: 


Tier  = 


(156) 
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with  amplifications: 


(157) 


A„cr(0)  =  eRn  and  Bn„( 0)  =  TjeRn/Rn 


where  Rn  =  r/n 


2  _  y^rjj 
cr  ~  4ct2/) 


and  tj,  the  final  circumferential  strain  is: 


Tf 


y/3pV02 

Aa 


(158) 


For  k  approaching  a  value  of  1/2  for  the  short  shell,  plane  stress  case,  the  strain  hardening 
moment  becomes  negligible,  giving  Sn  ~  0  and  i?2  ~  2 Tjn2.  The  maximum  harmonic,  nCT , 
is  determined  from  maximizing  equation  (152)  with  respect  to  n.  Reference  [4]  presents  the 
equation  for  nCT  resulting  from  a  graph  of  Qn  vs  Rn  as: 


3  _  c\/2t/(2  —  k)I\ i 
""  ""  3Fq2 


(159) 


where;  c=5/14  for  An  and  c=5/17  for  Bn.  Approximate  amplifications  for  the  critical  mode 
for  a2/?n2  <<  l,i?n2  ~  2 tjti2  and  k=  1/2  are  given  by  reference  [20]  as: 


_  gAn  __  T  f 

Ancr  =  Jjf  and  Bncr  =  (160) 

An  important  conclusion  from  studying  the  effect  of  length  in  plastic  flow  buckling  is  that 
the  number  of  waves  (n^)  in  the  buckled  mode  increases  with  shell  length.  For  the  same  load 
and  a/h  value,  a  short  shell  with  k=0.5  buckled  with  wave  number  13,  compared  to  its  longer 
counterpart  of  k=0,  which  buckled  with  wave  number  21  [21]. 

For  the  general  case  where  both  directional  and  strain  hardening  moments  are  significant, 
Lindberg  and  Florence  [4j  again  derive  an  expression  for  the  critical  mode  ncr,  by  minimizing 
An(£)  with  respect  to  n  to  obtain: 


4  -  t> 


(161) 


where:  C(q)  =  c(c xe  q,qo)  and  c=0.36,  cx  =  0.85,  <fo=0.1,  q  >  0.01,  q  =  2(2-t?K1T/fl»  /C1^)  ~ 
and  n'2  ^  (2 -%*0- 

The  approximations  to  the  amplification  functions,  A„(0),  can  be  used  to  establish  critical 
buckling  velocities  and  impulses  by  establishing  a  limiting  value  of  An(0).  This  is  done  in 
Section  3. 
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2.4.2  Elastic  Pulse  Buckling  of  Finite  Length  Shells 

Mclvor  and  Lovell  [23],  address  the  topic  of  autoparametric,  elastic  response  of  pulse 
loaded  finite  length  cylinders  in  the  same  manner  as  Goodier  and  Mclvor  [10]  for  infinite 
length  cylinders,  described  in  Section  2.2.  Expressions  for  the  membrane  and  bending  energy 
are  derived  and  given  in  Reference  [23].  These  are  considerably  more  complex  than  those  of 
equation  (42),  as  the  axial  strain  must  also  be  considered  in  the  energy  integral.  Displacement 
functions: 


u(x,0,  t)  =  y~]  Umn(t)  sin  COS  Ti  0 


m= 0  n=0 

oo  oo  _  _ 

/  \  \ — ^  r — r  .  ,  777  IT  X 

V{x,6,t)  =  Vmn(Q  COS— r— Sinng 

m=0  n=0 


tt’(x,0,<)  =  51  £W'mn(0cos^cos"tf 

_ _ _ n _ r\  * 


(162) 


m=0  n=0 


are  assumed  which  now  consider  axial  variation  of  the  shell  generator.  Boundary  conditions 
which  facilitate  a  solution  have  been  assumed;  §y,u,  radial  and  circumferential  shear  are  zero 
and  v  and  w  are  non-zero  at  the  ends.  These  boundary  conditions  are  not  neccessarily  of  most 
interest  in  practical  problems.  Expressions  (162)  are  substituted  into  the  energy  expressions 
which  are  in  turn  used  in  the  Lagrange  equation,  (47),  to  determine  the  equations  of  motion 
with  Umn.Vmn  and  Wmn  as  generalized  coordinates.  These  expressions  are  similar  in  form  to 
equations  (48)  and  (49)  of  the  infinite  length  case  but  are  considerably  more  complex  and  have 
three  pairs  of  equations  for  each  of  u,  v  and  w,  whereas  the  infinite  length  case  only  considered 
the  radial  displacement,  w. 

For  a  radial,  impulsive  velocity  on  the  shell,  the  equation  of  motion  (Lagrange's  equation) 
for  the  fundamental  hoop  mode  (n=m=0)  is  reduced  to  [23]: 


ti’oo  +  woo  =  0 


(163) 


With  initial  conditions  moo(0)  =  0  and  u’oo(O)  =  aV0/c ,  where  c  =  E / p(  1  -  v2)  =  Ei/p  and  the 
initial  velocity  Vo,  is  much  less  than  c  for  elastic  response,  the  solution  to  equation  (163)  is: 


,  ,  a  Vo  . 

woo(t)  =  - sin  r 

c 


(164) 


where  r  =  ct/a.  Perturbations  of  the  higher  order  harmonic  terms  are  used  to  investigate 
stability  with  initial  velocity  imperfections  given  by: 


Wmn  (0  )  =  7 


a  Vo 


mn  for  Tmn  ^  1 

C 


(165) 
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Figure  14:  Mathieu  Stability  Curve  for  Finite  Length  Shell  (from  reference  [23]) 

Considering  products  of  perturbations  to  be  negligible,  the  Lagrange  equations  for  generalized 
displacements,  Umn^Vmn  and  Wmn  become: 

U mn  +  C\Umn  +  C^Ymn  +  C3W  mn  =  0 

Fmn  +  C2  Umn  +  C4Vmn  +  C3Wmn  =  0 
Wmn  +  C3Umn  +  C5Vmn  +  Wmn(Ce  +  C4Woo/a)  =  0  (166) 

where  C,  are  defined  in  reference  [23]  and  are  functions  of  m,n,h,a  and  l. 

Neglecting  tangential  inertia  ( Umn  =  Vmn  =  0)  allows  Umn  and  Vmn  to  be  expressed  in  terms 
of  Wm „,  from  which  the  latter  of  equations  (166)  can  be  written  in  the  Mathieu  differential 
equation  form: 

limn  +  (ft  +  jU  sin  r)Wmn  =  0  (167) 

where  12  =  Ce,  -  %  +  and  is  a  function  of  the  shell  parameters,  m  and  n,  and 

^  =  C4Vo/c.  The  resulting  Mathieu  stability  curve  is  shown  in  Figure  14  for  the  ft  =  1/4  region 
of  stability  for  given  initial  velocity  and  shell  parameters.  Parametric  points  (ft,^)  are  defined 
by  intersections  of  families  of  curves  of  m  and  n.  Only  the  vn— 0  line  was  present  for  the  infinite 
shell  case  in  Figure  5.  As  a  result  of  finite  length  allowing  higher  axial  modes  of  response, 
there  are  many  more  points  in  the  unstable  region  of  the  Mathieu  diagram.  Reference  [23] 
demonstrates  that  as  the  shell  length  increases,  the  m=l,2,...  lines  converge  with  the  m=0  line, 
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giving  the  response  of  the  infinite  shell  (Figure  5).  The  response  of  the  shell  i  determined  by 
numerically  integrating  the  equations  of  motion  for  a  few  of  the  most  amplified  modes.  These 
modes  are  obtained  from  the  Mathieu  stability  diagram  and  confirmed  b  ensuring  that  the 
response  for  the  chosen  modes  is  of  comparable  energy  to  the  input  kinetic  energy.  Results  of 
numerical  integration  in  reference  [23],  show  that  the  modes  of  maximum  amplification  are  not 
always  of  the  m= 0  mode  as  in  infinite  shells. 

Since  this  analysis  is  for  purely  elastic  behaviour,  elastic  buckling,  which  grows  to  permanent 
deformation  or  collapse  before  hoop  mode  yield,  will  only  occur  for  very  thin  shells.  The 
important  aspect  of  this  theory,  as  in  the  case  of  the  infinite  shell  elastic  response  theory  is 
that  stresses  are  significantly  amplified  by  the  flexural  modes  over  the  purely  radial  hoop  mode 
stresses.  These  amplifications  need  to  be  determined  from  the  summation  of  the  harmonic 
amplitudes  obtained  from  numerical  integration  of  the  Lagrange  equations.  Results  reported 
in  reference  [23]  indicate  that  amplification  factors  as  high  as  four  for  circumferential  stresses 
and  eight  for  axial  stresses  may  exist. 

The  case  of  parametric  instability  of  elastic,  finite  length  cylinders  has  been  addressed  by- 
several  authors  for  nonimpulsive  loads  [19,15,16].  For  these  cases,  the  Donnell  shell  equations 
have  been  used  to  derive  the  equations  of  motion  for  the  unperturbed  radial  and  perturbed 
flexural  motions.  Bienick,  Fan  and  Lackman  [16],  derive  equations  of  motion  from  Donnell 
equations  and  use  the  Galerkin  method  to  obtain  a  solution  for  a  system  of  differential  equations 
and  arrive  at  a  Mathieu  type  equation.  Similar  results  to  those  of  Mclvor  and  Lovell  [23]  are 
attained,  where  an  increase  in  the  number  of  parametric  points  in  the  Mathieu  diagram  results 
from  decreasing  the  length  of  the  shell.  Yao  [15],  gives  a  good  discussion  of  the  development  of 
the  Mathieu  equation  from  Donnell  shell  theory  and  gives  several  examples  of  its  use  for  shells 
loaded  radially  and  axially.  It  is,  however,  for  very  thin  shells  under  static  or  periodic  loading. 

Anderson  and  Lindbeig  [19],  u'-e  Donnell  shell  theory  to  derive  the  equations  of  motion  for 
the  shell,  which  wall  be  repeated  here.  Using  the  compatability  equation: 


A  4F+  Eh[ 


Id2 


U’ 


d2 


a  dx2  '  dxdy 


w  l2  ,d2u\,d2w 

)  +fe)fe)l  =  ° 


dx2  A  Oy2 


(168) 


and  equilibrium  equation: 


1  d2 F  d2Fd2w  d2 F  d2w  d2F  d2u- 
a  dx 2  dy 2  dx2  +  ^  dxdy  dxdy  dx2  dy2 


(169) 


of  Donnell  shell  theory  [24],  inertial  and  initial  displacement  imperfection  terms  are  added,  from 
which  the  equilibrium  equation  becomes: 


>('”+"’4+2  ~£gS(w 


■(, 
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with  Nx  and  Nxg  defined  by  the  Airy  stress  function,  F,  and  Ng  defined  by  the  Airy  stress 
function  plus  the  membrane  force  from  radial  motion  giving: 


„  82F  „  d2F  ,  Eh  w0  .  d2F 

E  x  —  7aa7  >  ™ x9  —  rj/it,  and  A  $  —  .  ~  "h  o  2 

alo6 2  aodox  1  -  /nu2  a  ax2 


The  compatability  equation  becomes: 


A,F_Ehd^w 

a  dx2 


(171) 


(172) 


where  the  variations  of  w0  with  6  dissappear,  as  w  is  the  radial  hoop  mode.  The  stress  func¬ 
tion,  F,  represents  the  perturbed  flexual  motions.  Nondimensional  quantities  are  introduced  in 
reference  [19]  as: 

u  =  w/a,  =  W{/a,  £  =  x/a,  l  =  L/a,  and,  r  =  ct/a  (173) 

and  the  displacement  and  pressures  are  assumed  to  be  of  series  form: 

OO 

n(r)  cos  nd  sin  7 r£// 

n=l 

00 

*«.»)  =  E  <f>„(r)c os  n0sin  tt£/1 


P(f,0>r)  =  +  Z)p»(r)cos  nff  sin  tr£// 


(174) 


where  displacement  and  pressure  distribution  have  been  assumed  to  be  half  sine  wave  form 
(simply  supported  ends),  and  pn{T)  are  the  pressure  perturbations  from  which  stability  can  be 
determined.  Substituting  the  first  of  equations  (174)  into  equation  (172)  yields: 


1  ,  d2  d2  Eh  ^  7r2 

?( w +  w )F = — r  ^  7?”" tos nS sm 

’  n  =  l 


F  must  then  be  of  the  form: 


F  —  7 nt*n(r)  cos  n6  sin  7tf // 


(175) 


(176) 


which  when  substituted  into  equation  (170)  yields  the  equations  of  motion: 

fio  +  n0  =  po 

r  7/  2  7r\2  (1  —  /nii2)(7r//)4  2  ,  2  r  ,  . 

un  +  [a2(n2  +  —  )2  -(-  ^2  —  (-7r2//2)j:2"  “  n  uo\un  =  Pn  +  n2u0fn  (177) 

which  must  be  evaluated  numerically  to  determine  the  shell  motion.  Equation  (177)  is  of  the 
Mathieu  type. 
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2.5  Effect  of  Pulse  Duration  on  Response 


The  theory  of  the  preceeding  sections  has  been  for  response  to  ideal,  perturbed  axisymmet- 
ric  impulsive  loads  where  an  ideal  impulse,  applied  over  a  zero  time  period  can  be  represented 
as  a  Dirac  Delta  function.  The  ideal  impulse  was  applied  as  an  initial  velocity  condition  since 
it  had  no  time  dependence.  In  reality,  an  ideal  impulse  is  difficult  to  achieve,  as  the  pulse  will 
have  some  finite  time  duration. 

A  time  dependent  loading  function  must  be  incorporated  directly  into  the  equation  of  mo¬ 
tion.  This  was  seen  in  Section  2.4,  where  the  Donnell  shell  equation  was  used  to  investigate 
stability.  Anderson  and  Lindberg  [19],  investigated  the  effect  of  loading  rate  by  considering  the 
plastic  flow  equation  for  quasi-impulsive  loads  and  the  Donnell  equation,  (177),  for  quasi-static 
loading.  This  implies  that  dynamic  pulse  buckling  occurs  when  the  shell  is  in  the  plastic  flow 
regime,  which  is  true  for  all  but  very  thin  shells,  and  that  static  buckling  occurs  in  the  elastic 
regime.  Reference  [19]  defines  the  impulsive  loads  by  triangular  and  exponential  functions, 
giving  I  =  PT /2  for  the  triangular  shaped  pulse  and  I  =  PT  for  the  exponential  shaped  pulse, 
where  T  is  the  exponential  time  constant. 

The  plastic  flow  equation  of  motion  is  defined,  including  time  varying  pressure  and  variable 
material  properties,  as  [19]: 


u  -f 


a2Et  c^u 
E  ddA 


a 2Et  a g  d2u  a6  ae ,  d2 w, 

+  ('  E  +  E^ dB2  +  E^1  +  P  E^U'  +  dd2  ^ 


(178) 


where  the  initial  displacement  imperfection  and  the  displacement  and  pressure  functions  with 
perturbations  are  defined  as: 


OO 

u,  (6)  =  ^2  6n  cos  nO 

n=  1 

oo 

u(9,t)  =  «0(r)+  Un  COS  rid 

n= 1 
oo 

p(0,r)  =  p0(r)  +  ^p„cosn0  (179) 

n=l 

resulting  in  the  equations  of  motion: 

(7q 

uo  +  *^r(l  +  Wo)  —  Po 

unE{n2 -l)[~^n2  ~^}un=  +  ^(n2  -  l)*n  (180) 

with  initial  conditions  un(0)  =  un(0)  =  0.  Solutions  of  equation  (180)  are  obtained  by  numerical 
integration  as  are  the  results  for  the  elastic  model  frerr  equation  (177).  Both  equations  are  of 
Mathieu  form. 
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Figure  15:  Peak  Pressure  versus  Impulse  Curves  (from  reference  [19]) 

By  arbitrarily  assuming  a  critical  imperfection  amplification  of  1000  as  the  buckling  thresh¬ 
old,  Reference  [19]  derived  curves  from  the  two  theories  as  a  function  of  impulse  versus  pressure, 
shown  in  Figure  15.  The  two  curves  are  hyperbolic,  approaching  two  asymptotes  of  a  critical 
impulse  for  high  peak  pressure  and  critical  peak  pressure  for  high  impulse.  The  critical  impulse 
approaches  the  ideal  impulse  from  the  plastic  flow  pulse  buckling  theory  of  section  2.1  and  the 
critical  pressure  approaches  the  elastic,  static  buckling  pressure  for  infinite  time  duration.  The 
intermediate  portion  of  the  curve,  where  both  elastic  and  plastic  behaviour  are  significant,  is 
not  defined  by  either  theory.  Through  numerical  calculations  it  is  demonstrated  that  for  inter¬ 
mediate  a/h  ratios,  there  are  two  fundamental  growth  modes,  one  of  higher  harmonics  for  the 
plastic  flow  behaviour  and  one  of  lower  harmonics  for  the  elastic  behaviour. 

Lindberg  and  Anderson  [19],  investigate  the  effects  of  various  parameters  on  the  curve  of 
Figure  15.  The  pulse  temporal  variation,  of  exponential  or  triangular  shape,  showed  a  maximum 
shift  of  35  percent  in  the  curve.  Changing  the  a/h  ratio  had  an  expected  effect  of  changing 
the  relative  proportions  of  the  two  curves.  For  thick  shells  (small  a/h),  the  hyperbolic  curve  of 
the  plastic  flow  theory  was  dominant  with  very  little  or  no  elastic  curve  in  existance.  For  thin 
shells  (large  a/h),  the  opposite  behaviour  occurred,  with  the  elastic  hyperbola  dominating.  A 
change  of  length  of  the  cylinder  affected  only  the  elastic  branch,  as  the  plastic  flow  theory  used 
was  for  infinite  length  cylinders. 

Approximate  formulae  for  the  hyperbolic  curves  have  been  developed  in  reference  [19],  from 
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Figure  16:  Peak  Pressure  versus  Impulse  Curves  (from  reference  [191) 
which,  for  a  given  shell,  critical  load  parameters  can  be  established.  The  curves  are  defined  by: 

[(P/PA)-1][(I/IA)-1]  =  1  (181) 

where  for  the  plastic  flow  curve: 

3 

Pa  =  Pt  =  -trv(/i/o)  ,  and, 

Ia  =  It  =  {MlK)ll*a{poyy/2{hlafl2  (182) 

and  for  the  elastic  case: 

PA  =  PE  =  0.92  E(a/L)(h/a)5'2 

:  .id,  I A  —  IE  =  bpca(h/a )2  (183) 

where  K  is  the  a/ Et  slope  of  the  shell  material.  The  form  of  curve  derived  from  these  is  shown 
in  Figure  16. 

Characteristic  pulse  time  durations  are  also  derived,  defining  which  theory  best  describes 
the  shell  behaviour: 

T  >  Tt  =  2 (a/c)^=S  plastic  flow  theory 

T  <Te  —  5.5 (L/c)\/afh  elastic  theory 

Tg  <  T  <  Tj  intermediate  strain  reversal  theory  (184) 
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Nachbar  [25]  investigates  the  impulse  to  cause  failure  in  an  elastic/perfectly-plastic  infinite 
cylinder,  as  a  function  of  pressure  versus  time  duration,  for  a  rectangular  pulse  shape.  Failure 
is  assumed  to  occur  at  a  specified  displacement,  wj,  and  motion  termination  occurs  at  tj  when 
w  =  0. 

The  equation  of  motion  for  the  radial  hoop  mode  is  given  as: 

Ld2w  n  ah0  ,,oeS 

pah-^-r-  =  pR  -  - -  ( 1  So ) 

oti  1  —  w 


where  R  =  a(l  -  u;).  Upon  introducing  nondimensional  notation,  equationlSo  becomes  : 


d2w 
dr 2 


=  p(  1  -  u)  - 


(186) 


After  some  manupulation,  reference  [25]  presents  an  expression  for  the  critical  impulse  under 
Dirac  Delta  function  loading  as: 


4  = 


(187) 


For  elastic  perfectly  plastic  material  which  defines  cr(w)  =  u  for  u  <  uy  and  u  -  uy  for  u  >  uy, 
the  critical  impulse  for  a  Dirac  Delta  function  becomes: 


(188) 


where  wj  <  0.1a. 

For  the  case  where  the  duration  of  the  impulse  is  less  than  the  time  of  response  spent  in 
the  elastic  regime,  ty ,  the  impulse  to  cause  failure  is  given  as: 

I,  =  wyyJ A}  +  1  (189) 


where  A /  —  tj  —  ty 
given  as: 


1).  The  ratio  of  finite  pulse  to  Dirac  Delta  function  impulse  is 


lL  _  ;  _  111 

Ijt  sin  T /2 


(190) 


where  T  <  2arccos[A/^/Ay  +  1]. 


For  the  case  where  the  pulse  duration  is  greater  than  ty, 


^-  =  i  =  —TJa2,  +  1  (191) 

If,  wy  v 
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Figure  17:  Critical  Pulse  Parameters  for  Axisymmetrically  Loaded  Cylinders  (from  reference 
[25]) 


where  T  >  2arccos[Af,/A3  +  1]  and  pc  is  the  pulse  amplitude.  The  minimum  pulse  amplitude 

'  V  * 

which  will  cause  failure  defined  by  to/,  is: 


Po  =  toy[ 


i±^i, 

2  +  A}J 


occuring  over  a  pulse  of  duration: 


T  —  Tmax  —  A/v/2  -f  A^  +  —  + 


arcsin 


1  +  A  j 


(192) 


(193) 


From  this  theory,  curves  of  impulse  ratio  to  pulse  duration  time  were  produced  from  which 
critical  rectangular  pulse  parameters  can  be  derived.  These  curves  are  shown  in  Figure  17  as  a 
function  of  A/,  which  is  the  time  spent  in  the  plastic  regime. 

The  theory  of  Reference  [25],  has  been  based  solely  on  tracking  the  mean  radial  displacement 
of  the  hoop  mode,  to.  No  flexural  perturbations  are  considered.  It  is  assumed  that  flexural 
buckles  will  form  at  a  given  u>,  and  that  the  longer  the  time  spent  in  the  plastic  regime,  the 
more  likely  it  will  be  that  flexural  buckles  will  form. 
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2.6  Effect  of  Spacial  Pulse  Shape 

In  many  practical  problems,  and  certainly  in  the  case  of  a  submarine  subject  to  shock 
loading,  the  pulse  wave  will  not  be  axisymmetric  and,  in  fact,  will  also  not  be  applied  to  the 
entire  shell  instantaneously.  The  main  restriction  in  assuming  axisymmetry  in  the  preceeding 
theory,  is  that  the  circumferential  membrane  stress  is  independent  of  circumferential  location. 
In  asymmetric  loading  functions,  this  is  not  the  case.  For  thicker  shells,  which  buckle  in  plastic 
flow  at  higher  harmonics,  it  has  been  postulated  [27,19,4],  that  if  the  circumferential  stress  can 
be  assumed  to  be  constant  over  a  few  buckled  wave  lengths,  then  the  maximum  circumferential 
stress,  ag,  can  be  used  in  the  axisymmetric  pulse  buckling  equation.  This  requires  that  a  finite 
element  analysis  of  sufficient  accuracy  to  predict  the  peak  membrane  stress  of  the  shell,  ag,  as 
a  function  of  time,  be  undertaken.  These  values  are  then  used  in  the  numerical  integration  of 
the  axisymmetric  equations  for  perturbed  motion. 

In  thinner  shells,  or  for  shells  with  longer  duration  pulses,  the  buckled  wave  lengths  will  be 
longer  (smaller  number  of  harmonics,  n),  and  it  is  less  likely  that  the  circumferential  stress  will 
remain  constant  over  one  wavelength.  This  also  depends,  in  both  thick  and  thin  shells,  on  the 
degree  of  asymmetry  in  the  loading  function  which  must  be  at  least  smoothly  varying. 

An  impulse  loading  of  the  distribution  described  by: 

I  =  Ic  cos  0 ,  for,  -  7t/2  <  6  <  tt/2 

=  0,  for,  7r/2  <  9  <  3tt/2  (194) 

was  investigated  in  reference  [4],  Experimental  results  indicated  that  buckling  took  place  over 
about  120  degrees  of  the  shell.  Studies  utilizing  a$  results  from  finite  element  analysis,  com¬ 
paring  peak  impulse  values  to  give  the  same  deflection  amplitude,  indicate  that  a  greater  peak 
impulse  is  required  for  the  cosine  distribution.  This  varies  from  about  20  percent  greater  for 
the  thick  shell-plastic  flow  case  to  100  percent  greater  for  the  thin  shell-elastic  case.  It  was  also 
determined  [4]  that  for  the  cosine  pulse,  the  a/h  ratio  at  which  the  behaviour  is  purely  elastic  is 
shifted  much  higher.  For  the  axisymmetric  pulse,  a/h=288  and  for  the  cosine  pulse  it  increases 
to  480. 
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2.7  Strain  Rate  Reversal 

In  the  plastic  flow  theory  the  assumption  was  made  that  the  shell  was  in  a  continuous  state 
of  compression.  In  fact,  as  the  buckles  grow  with  time,  the  outer  and  inner  fibers  of  the  shell 
will  see  some  reversal  of  strain  rate,  with  tension  zones  eventually  occurring  in  the  convex  parts 
of  the  buckles.  The  effects  of  strain  rate  reversal  on  the  simple  theory  outlined  in  section  2.1 
were  investigated  by  Lindberg  and  Kennedy  [26]  through  a  numerical  finite  element  computer 
code  which  modelled  layers  through  the  shell  thickness  to  obtain  correct  strain  distribution  and 
more  complete  results.  Two  main  differences  from  the  simple  theory  were  noted  in  this  study. 

The  time  to  final  response  was  shorter  in  the  finite  element  code  results,  as  the  additional 
loss  of  energy  from  the  hoop  mode  to  the  flexural  mode  was  modelled.  In  the  simple  plastic 
flow  theory,  it  was  assumed  that  the  hoop  mode  absorbed  all  of  the  initial  kinetic  energy  before 
motion  ceased.  This  shorter  response  time  resulted  in  the  simple  theory  overestimating  the 
buckled  amplitudes  as  they  were  allowed  more  time  to  grow. 

The  occurrence  of  strain  rate  reversal  caused  the  buckles  to  unload  elastically.  This  greatly 
increased  the  stiffness  of  the  buckles  as  the  elastic  modulus  instead  of  the  strain  hardening 
modulus  was  in  effect.  The  motion  of  the  higher  harmonics  was  significantly  curtailed  as  the 
higher  harmonics  have  higher  curvatures  and  thus  experience  strain  rate  reversal  earlier.  The 
result  was  that  buckling  occurred  at  significantly  lower  harmonic  numbers  than  the  simple 
theory  predicted.  Figure  18  shows  the  results  for  harmonics  n=20  and  n=30,  comparing  the 
finite  element  computer  code  theory  to  the  simpler  plastic  flow  theory.  The  predominant  har¬ 
monic  of  the  simple  theory  was  n=30,  and  of  the  code  results  was  n=20.  The  simpler  theory 
significantly  overexaggerated  the  response  based  on  using  the  incorrect  n=30  harmonic  but 
reasonable  agreement  was  attained  for  the  n=20  harmonic  between  both  methods. 


52 


2.8  Material  Property  Effects 

For  impulsive  loading,  the  effect  of  high  strain  rate  on  material  properties  may  be  signifi¬ 
cant.  It  is  well  known  that  the  yield  stress  of  a  material  may  increase  under  high  strain  rates. 
This  effect,  although  acknowledged,  has  not  been  specifically  considered  in  available  references 
for  pulse  buckling.  Increased  yield  stress  will  affect  the  elastic  response  theory  (Section2.2)  if 
Young’s  Modulus  is  changed.  The  effect  of  increased  yield  stress  on  the  plastic  flow  theory  will 
be  to  increase  the  membrane  flow  stress  value  of  equation(3).  For  more  complete  theory,  where 
variable  material  properties  are  included  via  numerical  integration  of  the  equations  of  motion, 
a  strain  rate  sensitive  material  curve  could  be  used,  if  it  were  known. 

Another  strain  rate  effect  is  that  of  viscoelasticity  and  viscoplasticity,  where  the  modulus 
can  be  very  dependent  on  strain  rate.  For  response  to  pulse  loading,  the  strain  rate  continuously 
decreases  w'ith  the  motion  as  it  comes  to  rest.  This  means  that  a  continuously  varying  modulus 
must  be  considered  in  the  analysis.  Lindberg  and  Florence  [4]  have  addressed  pulse  buckling 
response  of  finite  length  cylindrical  shells  of  viscoplastic  material.  This  follows  the  theory  of 
Section  2.4  for  plastic  flow  buckling  of  finite  length  shells.  A  relationship  between  generalized 
stress  and  generalized  strain  rate  for  viscoplastic  material  is  used: 


7 c  =  cr/(To  -  1 


(195) 


where  7  is  the  viscosity  constant  from  material  tests  and  a0  is  the  generalized  stress  at  infinites¬ 
imally  small  strain  rate.  The  equation  of  motion  for  the  unperturbed  radial  motion  is: 


2(2  -  k)  ,  .  -(2  -  k)a0 

tt'o  +  r,  2  —  ~  77 

3  a‘p  A2np 

which  varies  from  the  plastic  case  by  the  term  ao^^'o  replacing  £/,tr0. 
tion  (196)  is  given  as: 


(196) 

The  solution  of  equa¬ 


t’d  =  [Vo  + 


3a 

2lK~2 


][1  _  e-2(2-*> 


<707f  3  a2p 

3a V  2(2-  k)a 07 


3a  t 

2^AV 


(197) 


The  equation  of  motion  for  the  perturbed  flexural  motions,  is  given,  in  the  nondimensional 
parameters  used  in  equation  (148),  as: 


u  - 


a2  3/c2  4yVo, d*ii  #2u 

(2  -  k)(K^  +  +  ~ 


(198) 


in  which  the  second  term  varies  from  the  non-viscoplastic  case  of  equation  (148).  By  making 
the  substitutions  of  equations  (127)  to  (130),  the  governing  equation  becomes: 

«„  -  (Pn  +  Qn/0"n  -  R2nun  =  0  (199) 
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where:  Pn  —  an<^  Qn  =  f^-£)A',  *  Wlt^  solution  in  the  same  form  as  the  non¬ 

viscoplastic  case: 

**»(£)  =  +  Bn(()bn  (200) 

where  A„(£)  and  Bn(£)  are  the  amplification  functions  for  the  displacement  and  velocity  per¬ 
turbations  and  are  now  defined  in  terms  of  Kummer  functions  and  given  in  reference  [4,  pg 
246]. 

Experimental  work  reported  in  Reference  [4]  demonstrated  that  for  fully  annealed  1015 
steel,  the  effects  of  increased  strain  rate  were  increased  the  yield  stress  and  lengthened  plastic 
portion  of  the  material  curve.  A  relationship  for  increase  in  yield  stress  with  strain  rate  is 
postulated  by  Symonds  and  Bodner[28]  as: 

°l°y  =  1  -f-  (c/D)1/p  (201) 

where  ay  is  the  static  yield  stress  and  p  and  D  are  constants  dependent  on  the  material  (5  and 
40.4  respectively,  for  mild  steel). 

The  key  material  parameter  in  plastic  flow  buckling  is  the  tangent  modulus.  If  the  tangent 
modulus  is  assumed  to  be  constant  (linear  plastic  regime),  closed  form  solutions  are  attainable. 
The  tangent  modulus  is  not  usually  linear,  in  which  case  numerical  integration  of  the  equations 
of  motion  is  required.  A  description  of  the  mateiial  curve  used  in  much  of  the  work  is: 

~  =  a/ E  =  c,  for,  0  <  e  <  ey 

•E't 

=  £ y  -f  K{(  -  ey),  for,  e  >  (y  (202) 

which  describes  curves  of  the  form  shown  in  Figure  19. 
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3  Approximate  Formulae  for  Critical  Impulse  and  Critical 
Modes  of  Dynamic  Pulse  Buckling 

From  the  preceding  theory,  formulae  for  the  critical  mode  of  response  and  for  approxima¬ 
tion  of  the  critical  impulse  to  produce  pulse  buckling  have  been  derived. 

The  critical  mode  number  (predominant  harmonic  of  response)  for  plastic  flow  buckling  of 
infinite  cylinders  or  rings  is  given  as  (equation  24)  [9]: 

ncr  -  ^(s2  +  1)  (203 ) 

which,  if  1  is  neglected  with  respect  to  s2,  becomes:  n „  =  v/6f./^.  Stuiver  [11]  derives  a 
critical  mode  equation  for  the  same  case  using  the  combined  elastic- plastic  theory  (equation  92): 

npl  =  0.85\/5j(|f  )1/2[1  +  (/?M£)2]i/<  (204) 

h  Lh  (y 

which  gives  comparable  values  to  equation  (203).  For  plastic  flow  buckling  of  short  shells  where 
directional  moments  dominate  the  strain  hardening  moments,  Vaughan  and  Florence  [20]  give 
the  critical  mode  (equation  135),  as: 


Tier 


(205) 


For  variable  length  shells,  Lindberg  and  Florence  [4]  give  the  expression  (equation  (161)): 


<  =  C(q) 


V^(2 -*•)*!  n„ 
3fc2a2  71  n 


(206) 


where  the  parameters  are  defined  in  section  2.4. 

For  viscoplastic  material  response,  where  the  viscoplastic  moments  dominate  the  directional 
moment,  "eference  [4]  gives  the  critical  mode  as: 


n 


3 

cr 


,  3a  \  /(2  -  k)p 

87a2  V  av? 


(207) 


where  the  parameters  are  defined  in  section  2.8.  All  of  the  above  formulae  are  for  thick  shells 
where  the  response  is  considered  to  be  entirely  plastic. 

For  response  where  behaviour  is  entirely  elastic,  equation  (54)  gives  the  critical  mode  for 
response  without  permanent  buckling  as  [10]: 


n 


CT 


1.316,/- 

fX 


(208) 
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Fot  larger  initial  impulses  where  buckling  can  occur,  the  empirical  relation  (equation  (65)) 
gives: 


n4  — 


0.25  0.13V2 


+ 


(209) 


Stuiver  [11],  gives  a  value  for  the  elastic  critical  mode  as: 


a  /Vo 
c 


Tier  =  2  — t/ — . 


(210) 


To  establish  the  critical  impulse  which  will  cause  dynamic  buckling,  a  limit  of  displacement 
amplitude  to  define  the  point  of  buckling  must  be  defined.  With  simplifying  assumptions  on 
material  properties  and  critical  mode  numbers,  equations  for  approximate  critical  impulse  have 
been  derived. 

For  plastic  flow  buckling  of  infinite  cylinders  or  rings,  Lindberg  and  Florence  [4]  derive  an 
expression  for  the  impulse  required  to  give  a  20  fold  increase  in  initial  shape  imperfection  as: 

Icr  =  y/3pca\l^(-)2  (211) 

y  il,  n 

Lindberg  and  Florence  also  give  threshold  impulse  values  for  variable  length  shells  failing  by 
plastic  flow  buckling,  by  taking  the  limit  as  a  shape  amplification  of  100,  ie.  An  -  100.  For 
long  shells,  with  k=0,  this  gives: 

Icr  =  3.07v/^a(~)2  (212) 

a 

and  for  short  shells  with  k=l/2,  this  gives: 

Icr  =  2.46v^g(-)3/2  (213) 

a 

for  q  >  0.5  where  q  =  ^or  var‘ab'e  length  shells,  the  impulse  functions: 
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=  (li 


k‘ ’  K\ 
4KTcJc 


V/4f(  #  v  W  -a  )3/2(ln  y/TeAn  )3/4 


for  0.01  <  q 
frr  0.01  >  q 


(214) 


can  be  used,  by  setting  the  amplification,  .4n,  equal  to  100,  and  where  the  variables  are  defined 
in  Section  2.4.  For  viscoplastic  material  response,  Lindberg  and  Florence  give  the  impulse 
amplification  relation: 


I  = 


4  2(2  -  k) 
3l  9  A'2 


a  5 


(215) 
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for  which  An  may  be  taken  as  100  to  determine  critical  impulses. 

For  buckling  in  the  elastic  regime,  where  the  flexural  buckles  exceed  yield  before  the  hoop 
mode,  a  value  of  p=4  is  used  in  equation  (71)  to  give  an  amplification  of  initial  shape  imper¬ 
fection  of  65.  From  the  relation  for  impulse  to  velocity  of  /  =  phVo,  the  critical  impulse  for  a 
condition  of  p=4  is  given  [4]  as: 

lev  =  1.15  pca(^f  (216) 

where  the  condition,  a/h  <  1.15/cv,  must  be  met  to  ensure  the  hoop  motion  does  not  exceed 
yield. 

The  curves  presented  in  Section  2.5  from  reference  [19]  can  also  be  used  to  establish  critical 
impulses,  particularly  if  the  loading  is  not  lor  an  ideal  impulse. 

To  cover  the  full  range  of  a/h  values,  Lindberg  and  Florence  [4]  suggest  using  the  plastic 
flow  and  elastic  buckling  threshold  equations,  and  then  using  whichever  gives  the  lowest  value. 
When  these  are  plotted  on  a  graph  of  a/h  versus  I& ,  Figure  20  results  (for  the  given  aluminum 
shell).  The  two  equations  defining  the  two  buckling  lines  are  given  as: 

— —  =  1.807 D(h/a)3/2  for  a/h  <  . 405/D 2 
pea 

=  1.15(/i/a)2  for  a/h  >  .405/D2  (217) 

The  results  of  these  various  formulae  have  been  investigated  via  a  computer  program  for 
shells  of  various  dimensions  and  properties.  Table  2  gives  the  shell  parameters  and  Table  3 
gives  the  results. 

Models  1  to  19  were  used  to  investigate  shells  with  increasing  a/h  ratios.  Columns  A,  B  and 
C  give  the  critical  mode  from  infinite  cylinder,  plastic  flow  theory.  The  critical  mode  increases 
with  a/h  and  the  theory  becomes  inapplicable  for  a/h  much  greater  than  40,  where  elastic 
behaviour  starts  to  have  some  influence.  Columns  D  and  E  give  the  critical  modes  for  the  short 
cylinder,  plastic  flow  theory.  The  greater  resistance  to  bending  resulting  from  the  directional 
moment  contribution  is  evident  in  the  lower  mode  numbers  for  these  cases.  Columns  F,  G  and 
H  give  the  critical  impulse  and  velocity  to  cause  plastic  flow  buckling.  As  the  shell  becomes 
thinner  (larger  a/h),  the  required  impulse  becomes  less,  as  would  be  expected.  Column  I  gives 
the  critical  mode  for  the  low  velocity,  elastic  vibration  response.  Columns  J  and  K  give  the 
critical  modes  for  the  larger  velocity,  elastic  buckling  case.  The  critical  modes  also  increase 
with  a/h  ratio.  Column  L  give.,  the  critical  impulse  and  velocity  to  cause  elastic,  dynamic  pulse 
buckling  and  as  in  the  plastic  flow  case,  as  a/h  increases,  the  required  impulse  decreases. 

The  effect  of  changing  the  material  yield  stress,  oy,  is  investigated  with  Models  3,  20  and 
21.  There  is  no  change  for  any  of  the  formulae  except  for  column  B  which  shows  a  decrease  in 
the  critical  mode  number  with  decreasing  yield  stress.  This  is  from  the  elastic-plastic  theory  of 
Stuiver  [11]  who  included  some  elastic  effect  in  the  determination  of  the  plastic  flow  buckling 
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critical  modes.  Formulae  for  column  A  and  C  would  also  be  affected  as  they  are  a  function  of 
the  flow  stress,  am,  which  is  a  function  of  the  yield  stress.  This  can  be  seen  in  models  31  to  33 
where  am  is  varied. 

Models  3  and  22  to  25  were  used  to  investigate  the  effect  of  varying  the  tangent  modulus 
in  the  strain  hardening  portion  of  the  material  curve.  Columns  A,  B  and  C  for  the  infinite 
shell,  plastic  flow  theory  show  a  dramatic  decrease  in  the  critical  mode  for  increasing  tangent 
modulus.  Columns  F  and  G  show  that  the  required  critical  impulse  increases  with  tangent 
modulus  for  the  infinite  shell,  plastic  flow  case,  as  would  be  expected  due  to  an  increase  in  the 
shell  stiffness.  Columns  D  and  E  for  the  critical  mode  and  column  H  for  the  critical  impulse 
show  no  change  with  varying  tangent  modulus  as  these  values  are  for  the  short  shell  case  where 
the  directional  moment  is  assumed  to  be  dominant  over  the  strain  hardening  moment. 

The  effect  of  varying  the  initial  velocity  was  investigated  with  Models  3  and  26  to  30.  There 
is  a  slight  variation  in  column  B  with  the  critical  mode  number  increasing  with  increasing  initial 
velocity.  This  is  from  equation  (204)  which  shows  a  variation  of  n„  as  the  square  root  of  F0. 
Column  E  shows  a  strong  dependency  of  on  Vo.  The  other  plastic  flow  formulae  show  no 
variation  with  initial  velocity.  The  critical  impulse  and  velocity  formulae,  are,  of  course,  not 
affected  by  Vo-  The  elastic  critical  modes  for  velocities  large  enough  to  cause  buckling  are 
influenced  by  the  initial  velocity.  Columns  J  and  K  show  that  the  critical  mode  number  varies 
proportionally  to  the  initial  velocity. 

Models  3  and  31  to  33  were  used  to  investigate  the  effect  of  varying  the  flow  stress  value, 
trm.  Columns  A  and  C  show  a  decrease  in  critical  mode  with  decreasing  flow  stress.  Column 
B,  from  equation  (204),  shows  no  change  in  critical  mode  with  flow  stress,  but,  the  flow  stress 
is  an  average  value  which  is  dependent  on  the  yield  stress  and  tangent  modulus,  which  are 
parameters  in  equation  (204)  so  some  effect  should  be  realized.  This  is  also  the  case  in  the 
critical  impulse  formulae  of  columns  F  and  G  w'here  the  tangent  modulus  is  used.  The  flow- 
stress  does  not  have  any  influence  on  the  elastic  theory. 
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Model 


1 

2 

3 

4 

5 

6 


Min)  I  ^y(psi) 


60,000 
pr\  nnn 

rJ 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 

60,000 


5 


70,000 

70,000 

70,000 

70,000 

70,000 

70,000 


,000 

,000 


70,000 

70,000 

70,000 

70,000 


,0 

,0 

,0 


65,000 

61,000 

60,000 


Efc(psi) 


1,000,000 

i  one  onn 

*  V-  Vi  ,V> 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 
1,000,000 


1,000,000 

1,000,000 


2,000,000 

3,000,000 

4,000,000 

10,000 


1,000,000 

1,000,000 

1,000,000 

1,000,000 

1,000,000 


1,000,000 

1,000,000 

1,000,000 


Vo(in/sec) 


6, 

6; 

6, 

6, 

6, 

6, 

6, 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 
6,000 


6,0 

6,0 


p(lb-sec2/in4) 


0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 

0.000787 


0.000787 

0.000787 


0.000787 

0.000787 

0.000787 

0.000787 


0.000787 

0.000787 

0.000787 

0.000787 

0.000787 


0.000787 

0.000787 

0.000787 


Table  2:  Shell  Parameters  Used  in  Table  3 


Model 


Plastic  Flow  Theory 


Hi  i ::  1  :  c 


A 

B 

N 

C 

cr 

D 

E 

■^cr 

j-  F 

-Vcr( in/s) 
G 

H 

Ncr 

'  I  J 

K 

-  j 

I  -V 

Acr  ¥cr 

L 

1 

6 

5 

6 

9 

8 

[  5.1-6465 

8.6- 10943 

5.8-7336 

4 

4 

3 

18 . 5- 23537 

2 

12 

10 

12 

13 

14 

2.5  3232 

4.3-5471 

4.1-5187  M  5 

7 

6 

9 . 3- 11768 

A 

l? 

16 

_  18 

15 

18 

1.7-2155 

2.9-  3648 

3.3-4236 ; 

;  7 

11 

10 

6 . 2- 7846 

4 

25 

21 

24 

18 

22 

1.3-1616 

2.2-2736 

2.9-3668 ' 

8 

14 

13 

4.6- 5884 

5 

32 

26 

30 

20 

25 

1.0-1293 

1 . 7-2189 

2.6-3281  : 

j  9 

18 

17 

3.7-4707 

6 

38 

32 

36 

22 

29 

.8-1077 

1.4-1824 

2.4-2995  j 

!  10 

21 

20 

3.1-3923 

7 

45 

37 

•4 

42 

24 

32 

. 7-923 

1.2-1563 

2.2-27  73  ! 

1  li 

25 

23 

2.6-3362 

8 

51 

43 

48 

26 

35 

.6-808 

1.1-1368 

2.0-2594  | 

li 

28 

27 

2 . 3-2942 

9 

58 

48 

54 

27 

38 

.55-718 

1.0-1216 

1.9-2445  ! 

'  12 

32 

30 

2 . 1-2615 

10 

64 

53 

60 

29 

40 

.5-647 

.86-1094 

1.8- 2320  j 

1  13 

35 

34 

1.9-2353 

11 

77 

64 

72 

31 

46 

.4-539 

. 7-912 

1.7-2118  1 

!  14 

42 

41 

1 . 5-1961 

12 

90 

75 

84 

34 

50 

. 36  462 

.6- 782 

1.5-1961  | 

!  15 

49 

47 

1 . 3-1681 

13 

103 

86 

96 

36 

55 

. 32-404 

.54-  684 

1.4-1834  j 

'16 

57 

54 

1.2-1471 

14 
,  _ 

116 

96 

108 

39 

60 

.28-359 

.48-608 

1.4-1729  : 

17 

64 

61 

1.0-1307 

15 

129 

107 

120 

41 

64 

.25  323 

.43-547 

1.3-1640  : 

18 

71 

68 

.93-1177 

16 

162 

134 

150 

46 

74 

.2,  259 

.34-438 

1.2-1467  j 

!  20 

89 

85 

.  74-942 

17 

194 

161 

180 

50 

84 

.17-216 

.29-365 

1.1-1340  ! 

22 

106 

102 

.62-785 

IS 

226 

i»a 

211 

54 

93 

.15-185 

.25-313 

.98-1240  j 

24 

124 

119 

.  53-6  73 

19 

259 

215 

241 

58  101 

^ . 13-162 

.22-274 

.91-1160  ! 

26 

142 

136 

46- SftR 

20 

19 

15 

18 

15 

18 

“  1~  7-2155 

2.9-3648 

3.3-4236 

7 

11 

10  5 . 2- 7B46 

21 

19 

13 

18 

15 

18 

1.7-2155 

2.9-3648 

3.3-4236 

7 

11 

10 

T2 

A  A 

13 

“12 

T5 

T3 — 

18" 

2.4-3048 

471-5159 

3.3-4236 

7 

11 

10 

U • 4“ 7040 

6 .2-7846 

23 

O  A 

11 

11 

10 

15 

18 

2.9  3721 

5.5-6318 

3.3-4236 

7 

11 

10 

6.2- 7846 

9 

11 

9 

15 

18 

3.4-4301 

5.7- 7296 

3.3-4236  : 

!  7 

11 

10 

6.2- 7846 

25 

194 
— — 

152 

180 

15 

18 

.17-216 

.29-  365 

3.3-4236 

:  7 

11 

10 

6 . 2-  7846 

26 

19 

1/ 

18 

15 

21 

1.7-2155 

"2.9-  3648 

3.3-4236 

7 

13 

12 

6.2-7846 

2  / 

19 

15 

18 

15 

17 

1.7-2155 

2.9-3648 

3.3-4236 

7 

10 

9 

6 . 2-7846 

28 

19 

15 

18 

15 

14 

1.7-2155 

2.9-3648 

3.3-4236 

7 

8 

7 

6 . 2-7846 

29 

19 

15 

18 

15 

10 

1.7-2155 

2.9-3648 

3.3-4236 

7 

7 

4 

6 . 2-  7846 

3  0 

19 

15 

18 

15 

8 

1.7-2155 

2.9-3648 

3 . 3-4236 

7 

7 

9 

f.  9  70  AC 

31 

18 

16 

17  '* 

15 

18  ' 

"l .  7-2155 

2.9-3648 

3.2-4082 

7  ~ 

Tf 

10 

D  .  £.  —  /  O  4 

~6 .2-18  46 

32 

18 

16 

16 

15 

18 

1. 7-2155 

2.9-3648 

3.1-3954 

7 

n 

10 

6 . 2- 7846 

33 

18 

16 

16 

15 

18 

1.7-2155 

2 . 9-3648 

3.08-3921 

7 

n 

10 

6.2- 7846 

A=  Eqn  20  3,  Re  f  9 

E=  Eqn 

206,  K=l/2,  Ref  4 

1  = 

Eqn 

208, 

Ref  10 

B=  Eqn  204,  Ref  11 

F=  Eqn 

211,  Ref 

4 

J= 

Eqn 

209, 

Ref  4 

C=  Eqn  206,  K= 

0,  Ref 

1  G=  Eqn 

212,  K=0. 

Ref  4 

K= 

Eqn 

210, 

Ref  11 

D=  Eqn  205.  Ref  20 

H=  Eqn 

213,  K=l/2,  Ref  4 

L= 

Eqn 

216, 

Ref  4 

Table  3:  Results  of  Various  Formulae  for  Critical  Modes  and  Impulses 
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4  Review  of  Numerical  Solutions  to  Dynamic  Pulse  Buckling 

Finite  difference  and  finite  element  numerical  methods  have  been  used  to  solve  a  wide  vari¬ 
ety  of  complex  structural  mechanics  problems,  particularly  where  analytical  solutions  cannot  be 
easily  found.  The  successful  utilization  of  numerical  methods  hinges  on  being  able  to  formulate 
a  discretized  model  and  solution  scheme  capable  of  modelling  the  correct  physical  behaviour. 
For  dynamic  pulse  buckling,  this  means  that  the  model  must  be  capable  of  reproducing  the 
nonlinear,  elasto-plastic  motion  of  the  shell  wall  in  the  higher  harmonic  modes  which  occur  in 
pulse  buckling. 

Most  numerical  studies  of  pulse  loading  consider  only  unperturbed  dynamic  response  and  do 
not  consider  the  growth  of  buckling  modes  in  the  solution.  Element  grid  discretization,  nonlinear 
formulations,  solution  methods,  material  property  representation  and  the  failure  cutoff  criteria 
are  parameters  which  need  to  be  investigated  for  numerical  determination  of  dynamic  buckling. 

Ishizaki  and  Bathe  [29],  investigated  static  and  dynamic,  linear  and  nonlinear  response  of 
perfect  and  imperfect  shells  with  the  finite  element  program  ADINA  [30].  A  spherical  cap.  a 
cylinder  and  a  sphere  were  investigated  to  determine  the  collapse  loads.  For  large  displacement, 
elasto-plastic  (E.P.)  response,  updated  and  total  Langrangian  (T.L.)  kinematic  formulations 
were  investigated  as  well  as  the  modified  Newton  and  the  BFGS  (Broyden-Fletcher-Goldfarb- 
Shanno)  iterative  solution  schemes  for  equilibrium  [31].  Without  modelling  imperfection  per¬ 
turbations  in  the  shell  geometry,  the  elastic  bifurcation  and  the  elasto-plastic  yield  collapse 
loads  can  be  determined.  Figure  21a  demonstrates  the  response  of  a  complete  sphere  to  a  static 
loading.  Figure  21b  shows  the  response  of  the  perfect  sphere  to  a  dynamic  step  loading  of  half 
the  elastic  buckling  k  magnitude.  The  elasto-plastic  curve  assumes  a  permanent  set  about 
which  the  response  oscillates  and  the  elastic  response  oscillates  with  greater  amplitude  but 
reaches  zero  displacement  on  each  oscillation.  Figure  21c  shows  the  response  to  static  loading 
of  an  imperfect  shell.  The  lower  curves  include  the  geometric  nonlinear  total  Langrangian  for¬ 
mulation  which  allows  the  growth  of  the  imperfections  to  produce  instability.  Figure  21d  shows 
the  response  of  the  imperfect  shell  to  the  dynamic  step  load.  The  top  curve  shows  the  results 
of  the  elasto-plastic,  total  Langrangian  formulation  which  models  the  unbounded  growth  of  the 
imperfection.  This  latter  curve  models  the  physical  characteristics  that  the  analytical  solutions 
did,  that  is,  allowing  unbounded  nonlinear  growth  of  displacements  for  a  specified  load.  Two 
important  factors  in  reproducing  this  behaviour  with  the  finite  element  method  are  that  im¬ 
perfections  in  the  modes  of  response  must  be  modelled  to  produce  the  unbounded  growth,  and 
that  the  load  amplitude  has  to  be  increased  in  consecutive  analyses  until  instability  occurs. 
Instability  was  established  in  the  ADINA  study  at  the  point  where  the  determinant  of  the 
stiffness  matrix  became  negative  (singular).  Some  major  difficulties  in  formulating  a  solution 
to  pulse  buckling  response  of  shells  with  the  finite  >  ’ement  method  are  evident. 

The  analytical  solutions  discussed  in  Section  2  indicate  that  the  circumferential  modes 
of  predominant  response  (buckling  growth)  occur  in  higher  harmonics  <han  static  buckling 
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PRESSURE  p/p 


modes.  It  can  also  not  be  easily  determined  which  mode  will  predominate.  Therefore  it  is 
necessary  to  model  imperfections  for  several  harmonics  of  numbers  up  to  30  or  higher.  This 
requires  a  very  high  level  of  discretization  in  the  circumferential  direction  of  the  shell.  In 
determining  the  response  of  the  cylindrical  shell  to  a  cosine  pulse  (the  same  problem  as  addressed 
by  Lindberg  and  Kennedy  [26]),  the  model  was  discretized  with  60,  8  noded  elements  in  one  half 
of  the  circumference.  Models  to  investigate  dynamic  pulse  buckling,  particularly  as  length  and 
complexities  such  as  stiffeners  are  included,  will  very  quickly  attain  a  large  number  of  degrees 
of  freedom.  To  determine  the  limit  load  of  dynamic  pulse  buckling,  several  analyses  runs  at 
increasing  load  levels  will  be  required.  Thus  determination  of  dynamic  pulse  buckling  response 
by  the  finite  element  method  will  be  an  expensive  and  time  consuming  proposition.  Quoting 
from  reference  [29],  ‘nonlinear  dynamic  buckling  analysis  is  frequently  beyond  the  current  state 
of  the  art’. 

Lindberg  and  Kennedy  [26],  investigate  plastic  flow  buckling  theory  with  a  finite  element 
code,  SABOR/DRASTIC  6  [32].  This  is  an  uncoupled,  axisymmetric  finite  element  analysis 
where  response  is  investigated  one  mode  at  a  time.  The  results  of  this  work  in  relation  to  the 
plastic  flow  theory  were  discussed  in  Section  2.7.  The  formulation  of  the  finite  element  problem 
involved  modelling  harmonic  imperfections  of  the  applied  impulse.  In  this  case,  harmonics 
derived  from  the  plastic  flow  theory  and  from  experimental  results  were  used.  Good  agreement 
between  the  numerical  results  and  experiment  was  attained,  due  in  part  to  the  fact  that  the 
imperfections  used  in  the  numerical  analysis  were  derived  from  the  experimental  results. 

Wesenberg  [18]  investigates  dynamic  buckling  response  for  shells  with  s^'.  s;al  a/h  ratios 
and  compares  them  to  experimental  results.  These  arc  A"-:.ussed  in  Section  2.3.  Here  again, 
imperfection  modes  and  amplitudes  were  measured  from  experimental  results  and  used  in  the 
numerical  finite  difference  formulation. 

More  recently  (1987-88),  Gefken,  Kirkpatrick  and  Holmes  [33,34,35]  applied  three  dimen¬ 
sional  nonlinear  finite  element  solutions  to  rings  and  finite  length  thin  shells  with  good  corre¬ 
lation  to  experimental  results.  The  necessary  requirements  of  using  a  finite  element  solution  in 
terms  of  initial  imperfect  shape  and  the  level  of  discretization  were  investigated  in  these  studies. 
It  was  determined  [35]  that  initial  imperfections  expressed  in  the  exponential  form: 

An  =  .05 h  n  <  10 

An  =  ~  n  >  10  (218) 

n1J 

gave  the  best  comparison  to  results  of  models  with  imperfections  derived  from  measuring  actual 
cylinder  imperfections.  The  finite  element  models  were  generated  with  a  series  of  imperfections 
covering  a  range  of  harmonics  such  that  the  geometry  was  formed  by  the  harmonic  summation: 

100 

R(6)  =  R  +  ^2  An  cos(n9  +  <j>n)  (219) 

71  =  2 
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Kirkpatrick  and  Holmes[34]  used  the  DYNA3D  finite  element  code  with  the  Hughes-Lui  shell 
element.  They  reported  requiring  approximately  ten  of  these  single-integration-point  elements 
per  buckling  wavelength. 

Many  studies  investigate  the  dynamic  response  of  shell  structures  to  impulsive  loading, 
but  do  not  investigate  the  possibility  of  dynamic  buckling.  Wu  and  Witmer  [36]  develop  a 
layered  finite  element  model  of  a  curved  beam  including  strain  hardening  and  strain  rate  effects. 
Comparisons  to  experimental  results  from  impulse  loading  are  made.  Lee  and  Horng  [2,37] 
develop  finite  difference  solutions  to  elasto-plastic  dynamic  response  of  ring  stiffened  cylinders 
to  shock  wave  type  loading.  The  critical  yield  points  occur  at  the  stiffener  shell  cor  nection. 
Buckling  could  have  been  investigated  if  the  harmonic  imperfections  were  included  in  tne  study. 

Several  studies  [38,39]  investigate  the  finite  element  solution  for  displacements  of  cylindrical 
shells  subject  to  underwater  shock  loads.  The  nonlinear  loading  function  of  the  coupled  fluid- 
structure  interface  is  modelled  using  doubly  asymptotic  approximation  and  boundary  elements. 
To  investigate  dynamic  stability  for  this  case,  a  nonlinear  finite  element  code  and  modelling  to 
include  imperfections  in  the  initial  shape  and  pressure  pulse  must  be  used  in  conjunction  with 
the  loading  algorithm. 


67 


5  Conclusions 


Analytical  solutions  for  specific  cases  of  dynamic  pulse  buckling  have  been  established  by 
various  authors  through  investigation  of  the  growth  of  perturbations  to  the  fundamental  motion. 
Cylinders  loaded  with  an  axisymmetric  pulse  which  respond  either  entirely  in  the  elastic  or 
entirely  in  the  plastic  material  range  have  been  investigated  analytically.  Approximate  formulae 
to  establish  the  critical  modes  of  buckling  and  the  critical  impulses  to  cause  buckling  for  various 
simple  cases  have  been  derived.  These  approximate  formulae  have  been  investigated  for  various 
shell  parameters  within  the  context  of  this  review.  The  dynamic  buckling  response  is  a  function 
of  the  shell  dimensions,  material  properties  and  loading  function.  No  analytical  solution  has 
been  derived  for  dynamic  buckling  of  shells  of  intermediate  a/h  ratios  which  encompass  many 
practical  cases. 

For  complex  shell  geometries,  complex  loading  functions  or  shells  of  intermediate  a/h  ratios, 
numerical  finite  element  or  finite  difference  methods  offer  a  potential  solution.  This  is  not 
straight  forward  and  few  studies  have  investigated  dynamic  buckling  solutions  using  numerical 
methods.  Material  and  geometric  nonlinearitics  have  to  be  included  in  the  formulation  and, 
as  in  the  analytical  solutions  for  the  simple  cases,  solution  is  for  excessive  growth  of  initial 
imperfections.  No  studies  of  ring-stiffened  shells  were  found  for  this  literature  review.  A  finite 
element  solution  seems  to  be  most  attractive  for  studying  dynamic  pulse  buckling  of  submarine 
structure.  Pressure  hulls  tend  to  be  of  an  intermediate  a/h  ratio  where  elasto-plastic  behaviour 
and  strain  rate  reversal  are  important  to  the  response. 

The  physical  concepts  of  dynamic  pulse  buckling  have  been  established  through  review  of 
the  analytical  studies.  These  and  the  approximate  formulae  for  critical  modes  and  loads  will 
be  invaluable  in  attempting  to  formulate  numerical  solutions  of  more  complex  problems  such 
as  pulse  buckling  of  ring  stiffened  submarine  pressure  hulls. 
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APPENDIX  A:  Derivation  of  the  Shell  Curvature  Expression 

The  change  in  curvature,  k,  is  defined  as  the  difference  between  the  deformed  shell  curvature, 
and  the  initial  shell  curvature, 

_  1  i_ 
p  a 

The  initial  length  of  element  mn,  from  Figure  A.l,  is  ds  =  ad<p,  and  the  initial  curvature 
is  The  length  after  deformation,  mjTii,  is  ds  +  Ads,  and  the  angle  of  curvature  is 

dcfi  +  A d<p,  which  gives  the  curvature  after  deformation  as: 

1  d(f>  4-  Ad(f> 
p  ds  +  Ads 

The  angle  A d<t>  is: 


Figure  A.l:  Geometry  of  Shell  Curvature  (from  reference  [40]) 
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The  change  in  length,  Adp,  is: 


ad<p  —  (a  —  w)  =  —i cdo  = 


ds 


giving  the  curvature  as: 

1  =  d<j>  +  {£$)ds 
P  ds(l-f) 

Neglecting  higher  order  terms,  this  reduces  to: 


1 

P 


1  ,  w 

-(!  +  -  + 
a  a 


d2w 


1  . ,  w %  1  02  w 

“(1  H - )  H — ?  IT77" 

u  a  a*  ap- 


which  when  substituted  into  the  expression  for  k,  gives: 


k  = 


1  d2w 
a2(  d<fi2 


+  w). 
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APPENDIX  B:  Derivation  of  the  Mathieu  Stability  Equation 

The  Mathieu  differential  equation  for  parametric  instability  can  be  most  easily  derived  for  the 
case  of  lateral  motion  of  an  axially  loaded  bar  [3].  The  differential  equation  for  the  lateral 
dynamic  motion  of  a  bar  under  axial  loading  is: 

d^w  r,d2w  d2w 

EIlh*  +  p  +  m~di T  =  0 

If  the  axial  loading  is  of  a  periodic  nature,  P(t )  =  PQ  +  P,  cos  <p/,  the  equation  of  motion  becomes: 

cP  w  d2  w  d2  u • 

EId^  +  ^Po  +  PtCOS^d^  +  m~&P=0 

If  the  response  is  assumed  to  be  periodic,  w(x,t )  =  fn(t)  sin  2Lp1.  the  equation  becomes: 


r  i  rr 

[m-^r  +  El 


4  4 

7*  7T4 


at2  '  p 

The  internal  expression  must  equal  0,  giving: 


fn  -  (Po  +  Pi  cos  <t>t) 


2  2 

n  n  ,  _  .  n~T 

fn]  sin  — —  =  0 


p 


1 


o2/  ,  2/i  po  +  Pt  cos  6t  x  ,  _  _ 

p ,  )fn  0  ,  72  —  1,2,3 

2  2  l  TP  J  .  1  1 

where  u>n  =  ~  is  the  free  vibration  frequency  of  an  unloaded  bar  and,  P*  =  ayj--£7  ,  is 

the  Euler  buckling  load  for  the  bar.  This  can  be  rewritten  as: 


fn  +  1  -  2 Hn  COS  Ot)fn  =  0 

which  is  the  traditional  form  of  the  Mathieu  equation  where  fin  =  sun \J^pfa ,  the  bar  frequency 

with  applied  axial  load  and  fin  =  2[P-’-p0)  •>  t^ic  excitation  parameter.  The  coefficient  of  fn 
approaches  zero  for  certajn  values  01  the  loading  function  resulting  in  instability.  The  more 
general  case  of  this  function  for  any  periodic  loading.  P(t)  =  P0  +  Pt$(t).  where  $(t  +  T)  =  <h(t ). 
is  known  as  the  Hill  equation: 


/  "  +  u2[\  -  2 /i$(0]/  =  o 

Functions  f(t),  must  be  found  to  satisfy  this  equation.  This  is  difficult  except  for  a  few 
simple  cases  (see  Bolotin[3]  and  McLachlin[14]  for  solutions).  Parametric  resonance  of  the  bar 
will  occur  if  the  excitation  frequency  is  twice  the  bar  frequency,  6  =  2fi.  Regions  of  instability 
can  be  plotted  as  functions  of  the  parameters,  J2  and  p.  The  shaded  regions  of  Figure  B.l  are 
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